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In  the  original  proposal  (Knight  1981),  the  general  goals  of  the  theo¬ 
retical  research  program  were  described.  Although  specific  programmatic 
changes  have  naturally  occurred  during  the  past  four  years,  the  fundamental 
objectives  remain  the  same,  namely  : 

o  To  determine  the  accuracy  of  a  theoretical  model  of  3-D  shock  wave- 
turbulent  boundary  layer  interactions  ("3-D  turbulent  interac¬ 
tions"),  where  the  theo^  ^tical  model  consists  of  the  3-D  mean  com¬ 
pressible  Navier-Stokes  equations  with  a  turbulent  eddy  viscosity 

o  To  investigate  the  flow  structure  of  3-D  turbulent  interactions  in 
simplified  geometries  through  a  close  cooperative  research  effort 
between  theory  and  experiment 

o  To  evaluate  the  hypothesized  physical  structure  of  the  3-D  turbulent 
interactions  at  a  variety  of  conditions  outside  the  range  of  experi¬ 
ments 

These  goals  represent  a  chronological  sequence  of  objectives.  A  major 
portion  of  the  first  three  years  focused  on  the  first  objective  (Knight  1982, 
1983,  1984a)  through  a  close  collaboration  with  the  Princeton  Gas  Dynamics 
Laboratory.  The  overall  evaluations  have  been  strongly  favorable,  and  provide 
the  impetus  for  achievement  of  the  second  goal.  During  the  fourth  year,  the 
research  program  has  made  significant  progress  in  achieving  the  second  objec¬ 
tive.  In  the  following  sections,  the  research  accomplishments  for  the  fourth 
year  and  research  program  for  the  fifth  year  are  presented. 
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Section  II.  Research  Accomplishments  for  the  Third  Year 
and  Research  Program  for  Fourth  Year 

A.  2-D  Turbulent  Interactions 

Although  the  principal  focus  of  the  research  effort  is  the  understanding 
of  3-D  turbulent  interactions,  a  modest  effort  has  been  directed  towards  2-D 
turbulent  interactions  over  the  past  four  years.  The  same  theoretical  model 
has  been  employed  for  both  the  2-D  and  3-D  research,  namely,  the  Reynolds- 
averaged  compressible  Navier-Stokes  equations  with  the  Baldwin-Lomax  turbulent 
eddy  viscosity  model. 

During  the  first  two  years  (Knight  1982,  1983),  the  objective  of  the  2-D 
turbulent  interaction  effort  was  the  examination  of  the  efficacy  of  the  theo¬ 
retical  model  through  comparison  of  the  computed  flowfields  with  the  mean 
flowfield  measurements  of  Settles  (Settles,  Gilbert  and  Bogdonoff  1980)  for 
the  2-D  supersonic  compression  ramp.  Calculations  were  performed  (see  Table 
1)  for  the  entire  experimental  matrix  of  Settles  at  a  nominal  Mach  number  of 
3,  ramp  angles  from  8  deg  to  24  deg,  and  Reynolds  numbers  ReSoo  =  0.76  x  106  to 
7.7  x  10^.  The  experimental  data,  consisting  of  surface  pressure,  skin  fric¬ 
tion,  surface  oil  flow  visualization,  and  boundary  layer  profiles  of  pitot  and 
static  pressure,  provided  an  extensive  database  for  examination  of  the  effi¬ 
cacy  of  the  theoretical  model.  The  numerical  algorithm  of  Beam  and  Warming 
(1978)  was  employed  for  solution  of  the  theoretical  equations.  The  con¬ 
clusions  of  the  research  effort,  sunmarized  in  Visbal  (1983)  and  Visbal  apd 
Knight  (1984),  are  the  following  : 

a.  The  determination  of  the  length  and  velocity  scales  of  the  Baldwin- 
Lomax  turbulent  eddy  viscosity  in  the  outer  region  is  unsuitable  in 
the  vicinity  of  separation.  The  model  predicts  a  sudden,  unphysical 
decrease  in  the  turbulence  length  scale  of  approximately  an  order  of 
magnitude,  resulting  in  an  corresponding  unphysical  reduction  in  the 
magnitude  of  the  eddy  viscosity. 


Table  1.  Theoretical  Research 
2-0  Shock  Wave  Turbulent  Boundary  Layer  Interactions 
Computations  of  2-D  Compression  Ramp 


Year 

Code  Mach  No . 
(nominal) 

Ramp  Angle  Reynolds  No. 

(ot,  deg)  (Reg00) 

Data  Comparison 

81-82 

B-W  3 

8  deg 

16  deg 

20  deg 

24  deg 

1.6 

x  106 

pw,  cf>  U,  M,  p 

pw,  cf,  U,  M,  p 

Pw*  cf-  xs-  xr •  u-  M 

pw,  cf,  xs,  xr,  U.  M 

82-83 

B-W  3 

20  deg 

0.76 

3.4 

5.6 

7.7 

X  106 

X  10® 

X  10® 

X  10® 

PW  ,  Xg  ,  Xj. 

P„.  xs-  xr 

Pw  xs*  xr 

Pw  xs>  xr 

B-W  2 

16  deg 

0.25 

x  10® 

'  xs  ’  xr 

83-84 

Development  of  2-D  Navier- 
hybrid  method  (MacCormack 

Stokes  code  using  MacCormack's 

1982) 

84-85 

MH  3 

8  deg 

16  deg 

20  deg 

1.6 

X  10® 

Pw,cf ,T,U,M,p 

Pw-Cf.T,U,M.p 

pw,cf ,T,xSIxr,U,M,p 

MH  2  16  deg  0.25  x  10®  p^  xs,  xr 


Note  :  There  are  typically  two  or  more  computations  for  each  case 

Legend  for  Numerical  Algorithm  : 

B-W  Beam  and  Warming's  Method  (Beam  and  Wanning  1978) 

MH  MacCormack's  Hybrid  Method  (MacCormack  1982) 

Legend  for  Data  Comparison  with  Experiment  : 

p^  Wall  static  pressure 

Cf  Wall  skin  friction  coefficient 

xs,  xr  Separation  and  Reattachment  points 
U  Velocity 

M  Mach  No.  Profiles 

p  Static  pressure  profiles 


a  c 


b. 


The  Baldwin-Lomax  model  exhibits  an  insufficient  upstream  propoga- 
tion  of  the  corner  interaction 


W 


c.  The  incorporation  of  a  ’relaxation  model ’  into  the  turbulent  eddy 
viscosity  model  improves  the  prediction  of  the  upstream  propogation. 
The  relaxation  model  introduces  an  additional  length  scale,  which  is 
observed  to  be  Reynolds  number  dependent. 

d.  The  Baldwin-Lomax  model,  with  or  without  the  relaxation  model,  fails 
to  accurately  predict  the  rapid  recovery  of  the  boundary  layer 
downstream  of  reattachment.  This  deficiency  was  attributed  to  the 
inability  of  the  theoretical  model  to  simulate  the  observed  rapid 
amplification  of  the  turbulent  fluctuations  (Settles,  Baca,  Williams 
and  Bogdonoff  1980,  Delery  1983)  across  a  shock-boundary  layer 
interaction.  The  effect  of  the  inclusion  of  the  relaxation  model  is 
to  diminish  the  turbulent  eddy  viscosity,  thereby  increasing  the 
upstream  propogation  as  noted  in  c.)  above.  While  improving  the 
flowfield  prediction  upstream  of  the  corner,  the  empirical  relaxa¬ 
tion  model  therefore  produces  the  wrong  behavior  downstream  of 
reattachment . 

On  the  basis  of  these  results,  a  second  effort  in  2-D  turbulent  interac¬ 
tions  was  initiated  during  the  second  year,  and  completed  in  the  fourth  year. 
The  objectives  of  this  effort  are: 

a.  Examination  of  the  sensitivity  of  the  computed  flows  to  the  numeri¬ 
cal  algorithm 

The  understanding  of  the  characteristics  of  a  numerical  algor¬ 
ithm  is  crucial  to  the  evaluation  of  the  theoretical  equations  which 
it  solves.  It  was  deemed  important,  therefore,  to  develop  a  second 
2-D  compressible  Navier-Stokes  code  using  the  popular  MacCormack 
hybrid  algorithm  (MacCormack  1982),  and  to  recompute  many  of  the 
same  cases  which  Visbal  had  computed  using  the  Beam-Warming  algor¬ 
ithm.  In  addition,  effort  was  focused  on  determining  whether  the 
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steady-state  numerical  solutions  obtained  using  the  MacCormack  hy¬ 
brid  algorithm  displayed  any  sensitivity  to  the  time  step  employed 
(i.e.,  Courant  number). 

b.  Direct  comparison  of  computed  and  measured  Reynolds  shear  stress  for 
the  2-D  compression  ramp 

Subsequent  to  the  theoretical  investigation  of  Visbal  and 
Knight,  a  series  of  experiments  (Muck  et  al  1983,  1984a, b)  were 
performed  at  the  Princeton  Gas  Dynamics  Laboratory  to  measure  the 
Reynolds  shear  stress  for  the  2-D  supersonic  compression  ramp  at 
Mach  3  for  ramp  angles  of  8,  16  and  20  deg  at  Regoo  =  1.6  x  10®. 
This  experimental  database  provided  the  opportunity  for  direct  com¬ 
parison  with  the  computed  Reynolds  shear  stress.  The  objective, 
therefore,  is  to  understand  the  reasons  for  the  failure  of  the 
Baldwin-Lomax  turbulent  eddy  viscosity  model  to  predict  the  rapid 
recovery  of  the  boundary  layer  downstream  of  reattachment. 

A  series  of  computations  were  performed  using  the  Reynolds-averaged 
compressible  Navier-Stokes  equations  and  the  same  Baldwin-Lomax  turbulence 
model  (see  Table  1).  Calculations  were  performed  for  a  Mach  2  compression 
corner  at  16  deg  and  Regoo  =  0.25  x  10®,  and  for  a  Mach  3  compression  corner  at 
8,  16  and  20  deg  and  Regao  =  1.6  x  10®.  The  computations  at  Mach  3  employed 
the  same  grid  spacing,  turbulence  model  parameters  and  upstream  profile  as 
employed  by  Visbal  and  Knight.  The  calculations  differed  only  in  the  numeri¬ 
cal  algorithm  utilized  (MacCormack  hybrid  vs.  Beam-Warming). 

The  conclusions  of  the  research  are  (see  Appendix  I)  : 

a.  The  computed  Reynolds  shear  stress  profiles  are  observed  to  be 
significantly  different  from  the  experimental  data.  The  height  of 
the  computed  peak  shear  stress  is  typically  a  factor  of  two  to  four 
too  small.  The  magnitude  of  the  computed  peak  shear  stress  is  in 
approximate  agreement  with  the  experimental  data  for  the  16  deg 
corner,  and  displays  a  modest  decrease  with  downstream  distance. 
For  the  20  deg  corner,  the  computed  peak  shear  stress  shows  a  rapid 
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decrease  with  downstream  distance,  in  disagreement  with  the  experi¬ 
mental  results. 


b.  The  Baldwin-Lomax  model,  based  upon  the  mixing  length  concept,  is 
incapable  of  accurately  predicting  the  recovery  of  a  separated  2-D 
compression  corner  flow.  The  mixing  length  model  is  formulated  on 
the  concept  of  an  equilibrium  turbulent  boundary  layer  exhibiting  a 
single  characteristic  velocity  scale  (Tennekes  and  Lumley  1972). 
Downstream  of  reattachment,  there  are  two  characteristic  velocity 
scales  of  the  turbulence,  namely,  1)  an  outer  velocity  scale  asso¬ 
ciated  with  the  turbulence  fluctuations  in  the  outer  portion  of  the 
reattaching  shear  layer,  and  2)  an  inner  velocity  scale  u*  = 
Aw(x)/pw(x)  associated  with  the  imposition  of  the  no-slip  boundary 
condition  downstream  of  reattachment,  which  creates  an  'inside  la¬ 
yer'  within  the  boundary  layer.  A  more  physically  realistic  turbu¬ 
lence  model  is  required  for  2-D  separated  compression  corner  flows, 
which  incorporates  the  effect  of  the  upstream  history  on  the  turbu¬ 
lence  and  the  oscillatory  motion  of  the  shock  wave  structure  (Dol¬ 
ling  and  Murphy  1982,  Dolling  and  Or  1983). 

c.  The  computed  flowfields  utilizing  the  Beam-Warming  and  MacCormack 
Hybrid  algorithms  are  overall  in  excellent  agreement.  This  implies, 
therefore,  that  the  effect  of  numerical  damping,  which  was  incor¬ 
porated  differently  in  the  two  methods,  is  negligible  on  the  flow- 
field  elements  examined  (e.g. ,  surface  pressure,  skin  friction,  and 
boundary  layer  profiles  of  velocity,  temperature,  density  and  Rey¬ 
nolds  shear  stress). 

• 

d.  The  steady-state  solutions  using  MacCormack’s  hybrid  algorithm  are 
observed  to  be  insensitive  to  the  Courant  number. 

It  is  evident,  therefore,  that  further  research  is  needed  to  provide  an 
adequate  qualitative  and  quantitative  understanding  of  the  separarated  2-D 
supersonic  compression  ramp  flows.  Firstly,  a  'new  look’  at  theoretical 
modelling  for  these  flows  is  required.  It  is  clearly  evident,  on  the  basis  of 
the  results  cited  above,  that  the  present  theoretical  model  (i.e.,  the  Rey- 


nolds-averaged  compressible  Navier-Stokes  equations  with  the  Baldwin-Lomax 
turbulent  eddy  viscosity  model)  does  not  embody  the  correct  physics  for  2-D 
separated  compression  corners.  Current  efforts  are  directed  towards  reviewing 
these  issues,  and  developing  alternative  approaches.  Secondly,  the  recent 
experimental  measurements  of  Reynolds  shear  stress  for  the  2-D  compression 
corner  by  Kuntz,  Amatucci  and  Addy  (1986)  using  an  LDV  differ  significantly 
from  the  earlier  measurements  of  Muck  et  al;  typically,  the  maximum  values  of 
the  kinematic  Reynolds  shear  stress  differ  by  a  factor  of  two  to  four  depend¬ 
ing  on  ramp  angle.  Additional  detailed  experimental  investigation  of  these 
flows  is  clearly  needed. 


B.  3-D  Turbulent  Interactions 

The  principal  focus  of  the  overall  research  program  is  the  understanding 
of  3-D  turbulent  interactions.  During  the  first  three  years  of  the  research 
effort,  the  principal  objective  was  the  determination  of  the  efficacy  of  the 
theoretical  model,  namely,  the  3-D  Reynolds-averaged  compressible  Navier- 
Stokes  equations  with  turbulence  incorporated  using  the  Baldwin-Lomax  alge¬ 
braic  turbulent  eddy  viscosity  model.  The  3-D  sharp  fin  configuration  (Fig. 
1)  at  Mach  3  was  selected  for  the  initial  investigations.  It  was  observed 
that  the  theoretical  model  provided  good  agreement  with  the  experimental  data 
(see  Table  2).  Consequently,  in  the  fourth  year  a  major  effort  focused  on  the 
development  of  the  mean  flowfield  model  for  the  3-D  sharp  fin  configuration. 
In  addition,  a  new  configuration,  the  3-D  swept  compression  corner,  was  com¬ 
puted  during  the  fourth  year.  In  the  following  sections,  the  results  of  the 
fourth  year  are  described. 

• 

1.  3-D  Sharp  Fin  at  20  deg 

The  configurations  computed  during  the  first  three  years  (Table  2)  in¬ 
volve  a  range  of  fin  angles  from  4  deg  to  20  deg  and  Reynolds  numbers  from  2.8 
x  105  to  9.3  x  10^.  Overall,  the  agreement  between  the  computed  and  experi¬ 
mental  results  is  very  good  (Knight  1984b,  1985a,  1985b;  Knight  et  al  1986) 
for  both  the  algebraic  turbulent  eddy  viscosity  of  Baldwin  and  Lomax  (Knight 
1984b,  1985a,  1985b;  Knight  et  al  1986)  and  the  two-equation  Jones-Launder 
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Table  2.  Theoretical  Research 


3-D  Shock  Wave  Turbulent  Boundary  Layer  Interactions 
Computations  of  3-D  Sharp  Fin 


Year 

Mach  No. 

Fin  Anale 

Reynolds  No. 

Data  Comparison 

(nominal ) 

(a,  deg) 

81-82 

3 

4 

9.3  x  105 

Pw* 

Pp. 

ch* 

sfc  visual 
yaw 

p,  pitch 

3 

10 

9.3  x  105 

Pw* 

pp* 

Ch, 

sfc  visual 
yaw 

p,  pitch 

82-83 

3 

10 

2.8  x  105 

Pw* 

Pp. 

sfc  visual 
yaw 

83-84 

3 

20 

9.3  x  105 

Pw* 

Pp. 

sfc  visual 
yaw 

► 


NOTE  :  Each  configuration  typically  represents  two  or  more  computations. 

The  purpose  of  multiple  calculations  is  to  investigate  the  effects 
of  grid  resolution. 


► 


► 


Legend  for  Data  Comparison  with  Experiment  : 


s?c  visual 


ch 

P 

pitch 


Wall  static  pressure 
Surface  flow  visualization  (oil  flow  or 
kerosene-lampblack ) 

Pitot  Pressure 
Yaw  Angle 

Wall  heat  transfer  coefficient  (Stanton  No.) 
Static  pressure  profiles 
Pitch  Angle 


► 
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(1972)  turbulence  model  (Knight  1985b,  Knight  et  al  1986),  for  which  the 
computations  were  performed  by  Dr.  C.  C.  Horstman  of  NASA  Ames.  On  the  basis 
of  this  verification  of  the  efficacy  of  the  theoretical  model  for  this  con¬ 
figuration,  the  computed  flowfields  using  both  theoretical  models  were  em¬ 
ployed  to  determine  the  structure  of  the  mean  flowfield.  The  principal  con¬ 
clusions  of  the  investigation  for  the  20  deg  sharp  fin  are  (Appendix  II)  : 

o  The  three-dimensional  velocity  fields  computed  by  both  theoretical 
models  are  in  close  agreement,  except  within  the  imnediate  vicinity 
of  the  surface  (specifically,  within  the  lower  approximate  20%  of 
the  boundary  layer),  where  differences  of  the  order  of  10%  to  20% 
exist  in  the  yaw  angle. 

o  The  computed  turbulent  Reynolds  shear  stress  profiles  differ  by  an 
order  of  magnitude  and  more  within  the  3-D  interaction  region 

o  The  similarity  of  the  computed  velocity  fields  and  significant 
difference  in  computed  turbulent  Reynolds  shear  stresses  imply  that 
the  3-D  sharp  fin  interaction  is  principally  an  inviscid,  rotational 
flow  except  within  the  immediate  vicinity  of  the  surface. 

o  The  computed  flowfields  display  a  prominent  vortical  structure  asso¬ 
ciated  with  the  shock-boundary  layer  interaction.  This  large  struc¬ 
ture  agrees  with  the  flowfield  models  of  Token  (1974)  and  Kubota  and 
Stollery  (1982).  The  calculated  flow  reveals  two  significant  sur¬ 
faces,  namely,  1)  a  three-dimensional  surface  of  separation  which 
originates  from  the  line  of  coalescence  (separation)  and  spirals 
into  the  center  of  the  vortical  structure,  and  2)  an  upper  surface, 
extending  upstream  into  the  undisturbed  portion  of  the  flow,  which 
defines  the  extent  of  the  flowfield  which  is  entrained  into  the 
vortical  structure. 

An  important  element  of  this  research  has  been  the  close  collaboration 
between  Dr.  C.  C.  Horstman  (NASA  Ames),  Prof.  S.  Bogdonoff  and  his  colleagues 
(Princeton  Gas  Dynamics  Laboratory),  and  the  present  author.  The  most  recent 
3-D  sharp  fin  configuration  (<*g  =  20  deg,  =  0.5  inch)  exemplifies  this 
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productive  interaction.  As  indicated  in  Section  III.B  and  III.C,  the  specific 
configuration  was  selected  after  detailed  discussion,  and  the  computations  for 
both  the  Baldwin-Lomax  model  (Knight)  and  Jones-Launder  model  (Horstman)  were 
performed  prior  to  the  experiment. 

2.  3-D  Swept  Compression  Comer 

The  second  major  focus  for  the  3-D  turbulent  interaction  research  in  the 
fourth  year  was  the  computation  of  a  selected  configuration  for  the  3-D  swept 
compression  corner.  This  specific  geometry,  described  in  terms  of  the  com¬ 
pression  corner  angle  <x  (measured  in  the  streamwise  direction)  and  sweepback 
angle  X  (Fig.  2),  represents  an  important  family  of  3-D  turbulent  interac¬ 
tions.  The  3-D  sharp  fin  and  2-D  compression  ramp  may  be  considered  specific 
cases  of  the  3-D  swept  compression  corner  family  corresponding  to  (a,x)  equal 
to  (90,<xg)  and  (<x,0)  deg,  respectively.  Extensive  experimental  data  has  been 
obtained  for  a  large  number  of  configurations  of  the  3-D  swept  compression 
corner  at  Mach  3  at  the  Princeton  Gas  Dynamics  Laboratory  (Settles,  Perkins 
and  Bogdonoff  1980;  Settles  and  Bogdonoff  1982;  McKenzie  1983;  Settles  and 
Teng  1984).  The  data  includes  surface  pressure  and  kerosene- lampblack  surface 
flow  visualization  for  more  than  forty  different  (<x,x)  configurations  at 
values  of  <*  and  x  up  to  24  deg  and  70  deg,  respectively  (Settles  1983). 
Detailed  flowfield  surveys  of  pitot  pressure  and  yaw  angle  have  been  obtained 
for  (<x,x)  =  (24,40)  deg  at  two  different  Reynolds  numbers. 

A  series  of  3-D  swept  compression  corner  configurations  have  been  calcu¬ 
lated  by  Horstman  (Settles,  McKenzie  and  Horstman  1984;  Horstman  1984)  using 
the  3-D  mean  compressible  Navier-Stokes  equations.  The  (<x,X)  =  (24,40)  deg 
configuration  was  computed  using  the  algebraic  eddy  viscosity  model  of  Cebeci- 
Smith  (1974)  and  the  two-equation  k  -  e  model  of  Jones-Launder  (1972)  for 
two  Reynolds  numbers.  The  computed  results  were  found  to  be  in  reasonable 
agreement  with  experiment,  although  discrepancies,  were  noted  in  the  surface 
pressure  distribution.  In  addition,  a  total  of  thirty-five  (35)  different 
configurations  were  calculated  by  Horstman  at  various  (ot.x)  using  the  Jones- 
Launder  (k-e)  turbulence  model  with  the  wall  function  model  of  Viegas  and 
Rubesin.  The  numerical  results  displayed  close  agreement  with  the  measured 
boundary  in  the  (a,x)  plane  between  cylindrical  and  conical  flow.  Recently, 
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however,  a  debate  has  arisen  concerning  the  existence  of  the  cylindrical- 
conical  boundary  (Wang  and  Bogdonoff  1986).  These  computations  displayed  good 
agreement  with  experimental  data  for  surface  pressure  for  a  number  of  (<x,x) 
configurations.  The  model,  however,  failed  to  accurately  predict  the  surface 
pressure  distribution  for  high  sweepback  angles  for  which  the  flow  exhibited  a 
large  separation  region  inferred  from  the  surface  streamlines.  In  particular, 
the  disagreement  for  the  (a,  X)  =  (24,  60)  deg  configuration  is  most  notice¬ 
able  (Fig.  3).  Subsequent  computations  by  Horstman  (unpublished)  utilizing 
finer  streamwise  grid  spacing  improved  the  prediction  of  the  surface  pressure 
within  and  downstream  of  the  corner  region;  however,  the  pressure  distribution 
upstream  of  the  corner  remained  essentially  unaffected. 


The  marginal  performance  of  the  Jones-Launder  turbulence  model  (using  the 
Viegas-Rubesin  wall  function  model)  for  the  (a,x)  =  (24,60)  deg  configuration 
motivated  a  reexamination  of  this  flowfield  using  the  theoretical  model  of  the 
present  author,  which  utilizes  the  Baldwin-Lomax  turbulence  model.  The  objec¬ 
tives  of  this  research  effort  are  the  following  : 


o  To  examine  the  accuracy  of  the  Baldwin-Lomax  model  for  the  class  of 
3-D  swept  compression  corner  interactions 


To  examine  the  sensitivity  of  the  computed  flowfield  to  the  turbu¬ 
lence  model  employed 


To  provide  a  theoretical  data  base  for  detailed  analysis  of  the 


flowfield  structure 


A  total  of  three  separate  computations  were  performed  for  the  (a,  X).: 
(24,  60)  deg  configuration.  In  all  cases,  the  freestream  Mach  number  M_  = 


2.95,  and  the  freestream  total  pressure  ptoo  =  100  psia. 


Case  1  :  Grid  System  and  Upstream  Boundary  Layer  Profile 
Equivalent  to  Horstman 


The  first  computation  of  the  (a,  x)  =  (24,  60)  deg  configuration 
employed  approximately  the  same  grid  spacing  as  Horstman  (1984).  A  total 


of  forty  (40)  streamwise  grid  planes  were  incorporated.  Within  each 
plane,  a  total  of  729  ordinary  grid  points  were  utilized  (27  x  27),  with 
ar.  additional  162  points  within  the  computational  sublayer.  The  total 
number  of  grid  points  was  35,640.  The  spanwise  grid  spacing  (z-direc- 
tion)  was  uniform.  The  grid  in  the  y-direction  was  stretched  geometric¬ 
ally  to  resolve  the  boundary  layer  on  the  surface. 

The  upstream  boundary  layer  profile  was  the  same  as  employed  by 
Horstman.  It  is  important  to  note  that  the  upstream  boundary  of  the 
computational  domain  was  parallel  to  the  corner  line  (Fig.  4)  for  both 
Horstman  and  the  present  author.  Horstman  utilized  the  same  boundary 
layer  profile  at  all  points  on  the  upstream  boundary  with  5^  =  0.38  cm), 
effectively  simulating  a  boundary  layer  which  had  developed  from  a  swept 
leading  edge.  In  the  experiment,  however,  the  leading  edge  of  the  flat 
plate  (to  which  the  swept  compression  corner  was  affixed),  was  perpen¬ 
dicular  to  the  oncoming  flow,  thereby  implying  a  boundary  layer  whose 
thickness  varied  on  the  upstream  boundary  of  the  computational  domain. 
This  effect  was  examined  later  (Case  2). 

The  freestream  total  temperature  Ttoo  was  495  deg  R,  and  the  wall 
temperature  520  deg  R  in  agreement  with  Horstman.  These  values  are  in 
agreement  with  the  experiment. 

Case  2  :  Modified  Upstream  Boundary  Layer  Profile 

The  second  computation  utilized  the  same  grid  system  as  Case  1,  with 
a  total  of  35,640  grid  points.  The  upstream  profile,  however,  was  chosen 
to  accurately  represent  the  oncoming  boundary  layer  in  the  experiment.  .A 
separate  boundary  layer  code,  developed  through  AFOSR  support  (Knight 
1983,  1984a),  was  utilized  to  compute  an  upstream  profile  corresponding 
to  the  flat  plate  leading  edge  perpendicular  to  the  oncoming  flow.  The 
boundary  layer  profile  on  the  upstream  surface  was  in  good  agreement  with 
the  experimental  boundary  layer  which  developed  on  the  flat  plate  in  the 
wind  tunnel.  The  boundary  layer  thickness  5^  on  the  upstream  boundary 
varied  from  0.091  cm  to  0.38  cm  over  the  13.2  cm  width  of  the  computa¬ 
tional  domain. 
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The  freestream  total  temperature  Ttoo  was  479  deg  R,  in  agreement 
with  the  experiment  (Kimmel  1985).  The  wall  temperature  505  deg  R  in 
agreement  with  the  experiment  and  previous  computational  studies  of  the 
3-D  sharp  fin. 

Case  3  :  Modified  Spanwise  (Z-Direction)  Grid  Spacing 

The  third  computation  utilized  the  same  grid  spacing  in  the  x-  and 
y-directions  as  Cases  1  and  2.  The  spanwise  (z-direction)  mesh,  however, 
was  modified.  As  indicated  previously,  the  spanwise  grid  spacing  for 
Cases  1  and  2  was  uniform  and  equal  to  0.508  cm.  At  the  apex  of  the 
swept  compression  corner  (i.e.,  at  the  intersection  of  the  corner  line 
and  the  plane  of  symmetry),  the  size  of  the  3-D  interaction  is  small.  It 
was  deemed  necessary,  therefore,  to  examine  the  effect  of  refining  the 
grid  in  the  spanwise  direction  in  the  vicinity  of  the  apex  in  order  to 
accurately  resolve  the  'inception  region'.  Consequently,  in  Case  3  a 
highly  stretched  grid  spacing  was  employed  in  the  z-direction,  with  a 
fine  spacing  az  =  4.57  x  10-3  cm  near  the  symmetry  plane  (z  =  0),  and  a 
larger  spacing  az  =  0.539  cm  near  the  right  boundary.  The  total  number 
of  grid  points  in  the  z-direction  increased  from  27  (Cases  1  and  2)  to 
46.  The  overall  total  number  of  grid  points  was  60,720. 

The  upstream  profile  was  identical  to  Case  2,  and  the  boundary  layer 
thickness  8^  on  the  upstream  boundary  varied  from  0.091  cm  to  0.38  cm 
over  the  13.2  cm  width  of  the  computational  domain. 

The  total  temperature  and  surface  temperature  were  identical  to  Ca?e 

2. 


The  results  for  Cases  1  to  3  are  presented  in  Figs.  5  through  7.  In  Fig. 
5a  to  5e,  the  computed  surface  pressure  is  displayed  at  five  spanwise  loca¬ 
tions,  corresponding  to  z  =  2.03  cm,  4.06  cm,  6.09  cm,  8.12  cm  and  10.15  cm. 
As  indicated  above,  the  boundary  layer  thickness  on  the  upstream  surface  (Fig. 
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4)  is  $3,  =  0.38  cam  for  Case  1  (in  agreement  with  the  profile  of  Horstman),  and 
varies  from  0.091  cm  to  0.38  cm  for  Cases  2  and  3.  The  vertical  axis  is  the 
surface  pressure,  normalized  by  the  upstream  static  pressure  p^.  The  horizon¬ 
tal  axis  is  the  streamwise  distance  x  -  xcorner.  where  xcorner  is  the  stream- 
wise  location  of  the  corner  line  at  the  specified  spanwise  location.  In  Fig. 
6,  the  computed  and  experimental  surface  pressure  are  displayed  in  conical 
coordinates,  where  the  horizontal  axis  is  (x  -  ^corner^2  +  zorigin^’  and 
zorigin  =  01,1  *s  the  z-coordinate  of  the  approximate  virtual  origin  as 

estimated  from  the  experimental  surface  flow  visualization.  The  experimental 
data  was  obtained  at  z  =  7.62  cm,  9.4  cm,  9.65  cm  and  10.16  cm.  The  computed 
profiles  in  Fig.  6  are  obtained  at  z  =  12.2  cm  Additional  profiles  (not 
shown)  indicate  that  the  computed  surface  pressure  is  approximately  conical 
outside  of  an  initial  inception  region  at  the  apex. 

Several  features  are  evident  from  Figs.  5  and  6  : 

1.  The  calculated  upstream  propogation  is  insensitive  to  the  boundary 
layer  on  the  upstream  surface. 

The  computed  surface  pressure  profiles  for  Case  1,  which  em¬ 
ployed  a  uniform  boundary  layer  thickness  on  the  upstream  boundary 
surface  (thereby  simulating  a  swept  leading  edge  to  the  flat  plate), 
is  quantitatively  very  similar  to  the  profiles  of  Cases  2  and  3, 
which  employed  a  non-uniform  boundary  layer  thickness  on  the  up¬ 
stream  boundary  surface  (thereby  simulating  a  straight  leading  edge 
to  the  flat  plate,  as  employed  in  the  experiment).  There  is  a  slight 
decrease  in  the  upstream  propogation  for  Cases  2  and  3  as  compared 
with  Case  1. 

2.  The  calculated  upstream  propogation  is  in  good  agreement  with  the 
experimental  data  as  indicated  in  Fig.  6. 

3.  The  computed  pressure  profile  upstream  of  the  corner  is  quantita¬ 
tively  different  from  the  experiment  (Fig.  6).  The  computed  profile 
fails  to  reproduce  the  rapid  pressure  rise  (although  the  location  of 
the  beginning  of  the  pressure  rise  is  accurately  predicted  as  in- 
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dicated  in  1.)  above).  In  addition,  the  drop  in  pressure  immediate- 
ly  upstream  of  the  corner  is  not  observed  in  the  computations. 

4.  The  peak  pressure  at  the  corner  is  accurately  predicted  (Fig.  6). 

5.  The  recovery  of  the  surface  pressure  downstream  of  the  corner  is 
somewhat  sensitive  to  the  upstream  boundary  layer  profile,  with 
Cases  2  and  3  displaying  a  more  rapid  recovery  than  Case  1.  Never¬ 
theless,  the  computed  pressure  downstream  of  the  corner  is  too  low. 

6.  The  computed  surface  pressure  is  insensitive  to  the  refinement  in 
the  grid  spacing  in  the  spanwise  direction  near  z  =  0,  i.e.,  the 
calculated  profiles  for  Cases  2  and  3  are  nearly  identical. 

In  Fig.  7a,  calculated  surface  pressure  profiles  are  displayed  in  conical 
coordinates  (x  -  xcorner)/(z  +  zorigin)  for  the  Baldwin-Lomax  model  (Cases  1 
to  3).  The  spanwise  location  of  the  profiles  is  typically  11.9  cm,  which  is 
within  the  region  where  the  computed  surface  pressure  profiles  display  approx¬ 
imate  conical  similarity.  In  Fig.  7b,  calculated  surface  pressure  profiles 
are  shown  for  three  cases  using  the  Jones-Launder  model  (Cases  4  through  6)  at 
a  spanwise  location  of  11.9  cm.  Case  4  represents  the  original  computation  by 
Horstman,  with  the  minimum  streamwise  grid  spacing  ax  =  0.38  cm.  The  computa¬ 
tions  using  the  Baldwin-Lomax  model  utilized  the  same  streamwise  grid  spacing. 
In  Case  5,  the  streamwise  grid  spacing  was  refined  with  ax^  =  0.25  cm.  Both 
Case  4  and  5  utilized  the  original  3-D  version  of  the  viegas-Rubesin  wall 
function  model  ('Viegas-Rubesin  I').  In  Case  6,  a  modification  of  the  Viegas- 
Rubesin  model  was  employed  ('Viegas-Rubesin  II').  In  Fig.  7c,  the  profiles  of 
Figs.  7a  and  7b  are  incorporated  in  a  single  plot.  Several  observations  are 
evident  from  Figs.  7a-c,  namely  .• 

1.  The  computed  upstream  influence  is  relatively  insensitive  to  the 
turbulence  model. 

2.  The  calculated  profiles  upstream  of  the  corner  consistently  fail  to 
reproduce  the  rapid  pressure  rise  and  subsequent  drop  (see  Fig.  6). 


3.  The  peak  pressure  at  the  corner  for  the  Jones-Launder  model  is 
sensitive  to  the  streamwise  grid  spacing  and  version  of  the  wall 
function  model  (i.e.,  Viegas-Rubesin  I  or  II). 

4.  The  peak  pressure  at  the  corner  for  the  Baldwin-Lomax  model  is 
insensitive  to  the  streamwise  grid  spacing.  The  computed  peak 
corner  pressures  for  Cases  1  to  3  and  6  are  essentially  identical. 

5.  The  pressure  profiles  downstream  of  the  corner  for  the  Baldwin-Lomax 
model  (Cases  1  to  3)  and  the  Jones-Launder  model  with  the  refined 
streamwise  grid  (Cases  5  and  6)  are  very  similar. 

At  the  present  time,  further  investigation  of  the  computed  flowfields  is 
in  progress.  Several  possible  causes  for  the  discrepancy  between  the  computed 
and  measured  surface  pressure  are  under  investigation,  including  a)  effect  of 
grid  resolution,  and  b)  inadequacies  in  the  turbulence  modelling. 

C.  Research  Program  for  Fifth  Year 

The  research  program  for  the  remainder  of  the  fifth  year  focuses  on  two 
principal  objectives  (Knight  1985c)  : 

1.  Analysis  of  3-0  swept  compression  Corner  : 

(a.  A)  =  (24,  60)  deg  at  Re^  =  2.5  x  105 

The  computed  results  for  the  3-D  swept  compression  corner  at  (a,  x)  = 
(24,60)  deg  and  ReSoo  =  2.5  x  10^,  obtained  during  the  previous  year,  are 
analyzed  in  detail.  A  detailed  comparison  of  the  computed  flowfields  using 
the  Baldwin-Lomax  model  (Knight)  and  Jones-Launder  model  (Horstman)  is  per¬ 
formed,  including  profiles  of  yaw  angle,  pitch  angle,  velocity,  Mach  number 
and  static  pressure.  A  comparison  of  particle  pathlines  is  also  performed  to 
examine  the  predicted  flowfield  structure.  A  significant  effort  is  focused  on 
the  examining  the  quality  of  the  numerical  grid  to  assure  the  fidelity  of  the 
computed  solution. 


2.  Calculation  of  3-D  Swept  Compression  Corner  : 

(a.  X)  =  (24.  6»)  deg  at  Re^  =  9.5  x  105 

The  3-D  swept  compression  corner  at  (a,  \)  =  (24,  60)  deg  and  Regoo  =  9.5 
x  10^  is  computed.  The  choice  of  this  configuration  is  based  upon  several 
factors  : 


a.  Concurrence  with  Experimental  Effort  at  Princeton  Gas  Dynamics 
Laboratory 

A  major  investigation  of  the  3-D  swept  compression  corner  at  (<*,  x) 
=  (24,60)  deg  and  Regoo  -  9.5  x  10-'  is  planned  for  1985-1936  at  the 
Princeton  Gas  Dynamics  Laboratory  (Bogdonoff,  Andreopoulos  and  Smits 
1985,  p.  66).  This  effort  will  include  a  detailed  flowfield  study 
near  feature  lines  (e.g.,  lines  of  coalescence)  and  in  the  incep¬ 
tion  region.  The  computation  of  this  configuration  during  the  same 
period  will  provide  opportunities  for  continued  close  interaction 
between  the  theoretical  and  experimental  efforts.  This  interaction 
can  include  the  utilization  of  the  computed  flowfields  to  suggest 
locations  for  experimental  measurements. 


b.  Complement  the  Previous  Calculation  of  the  3-D  Swept  Compression 
Corner  at  (a,  x)  =  (24,  60)  deg  for  Regoo  =  2.5  x  105 

The  present  calculation  will  complement  the  previous  computational 
study  of  the  same  geometry  at  the  lower  Reynolds  number  Regoo  =  2.5  x 
10^.  These  combined  studies  will  provide  a  detailed  examination  of 
the  effects  of  Reynolds  number  on  this  complex  interaction.  Prev¬ 
ious  experience  with  the  3-D  sharp  fin  at  ag  =  10  deg  indicated  that 
certain  features  (e.g.,  the  overshoot  in  pitot  pressure  upstream  of 
the  theoretical  inviscid  shock  and  outside  the  boundary  layer)  were 
accurately  predicted  at  the  lower  Reynolds  number,  but  not  as  close¬ 
ly  at  the  higher  Reynolds  number.  In  addition,  certain  features 
(e.g.,  the  overshoot  in  yaw  angle  in  the  same  physical  location) 
were  not  accurately  predicted  at  either  Reynolds  number.  This 


experience  strongly  suggests  the  importance  of  performing  separate 
computations  at  different  Reynolds  numbers  for  each  geometrical 
configuration  (i.e.,  specific  values  of  a  and  X)  of  the  3-D  swept 
compression  corner. 


iwinui'in  a  it!  w  mu  tutu  wutu  m  vu  ir’v1  ■  w"* 


Section  III.  Publications  and  Scientific  Interactions 


* 


» 


t 

I# 


A.  Written  Publications  -  Cumulative  Chronological  List 
1.  1  October  1981  -  30  September  1982 

a.  Knight,  D.,  "Application  of  Curvilinear  Coordinate  Generation 

Techniques  to  the  Computation  of  Internal  Flows",  in  Numerical 
Grid  Generation  -  Proceedings  of  a  Symposium  on  the  Numerical 
Generation  of  Curvi 1  inear  Coordinates  and  their  Use  in  the 
Numerical  Solution  of  Partial  Differential  Equations.  North- 
Holland,  New  York,  1982,  pp.  357-384.  [*]  [**] 

b.  Knight,  D.,  "A  Hybrid  Explicit- Implicit  Numerical  Algorithm  for 

the  Three-Dimensional  Compressible  Navier-Stokes  Equations", 
AIAA  Paper  No.  83-0223,  AIAA  21st  Aerospace  Sciences  Meeting, 
January  10-13,  1983.  Published  in  AIAA  J. .  Vol.  22,  Aug  1984, 
pp.  1056-1063.  [*]  [**] 

c.  Visbal,  M.,  and  Knight,  D.,  "Generation  of  Orthogonal  and 
Nearly  Orthogonal  Coordinates  with  Grid  Control  Near  Boun¬ 
daries",  AIAA  J. .  Vol.  20,  No.  3,  March  1982,  pp.  305-206.  [**] 
[***] 
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[*]  Research  sponsored  by  AFOSR  Grant  82-0040 

[**]  Research  sponsored  byAFOSR  Grant  80-0072 

[***]  Research  sponsored  by  AF  Contract  F-33615-C-3008 
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1  October  1982  -  30  September  1983 
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a.  Knight,  D. ,  "Calculation  of  a  Simulated  3-D  High  Speed  Inlet 

Using  the  Navier-Stokes  Equations",  AIAA  Paper  No.  83-1165, 
AIAA/SAE/ASME  19th  Joint  Propulsion  Conference,  Seattle,  Wash¬ 
ington,  June  27-29,  1983.  [*] 

b.  Visbal,  M. ,  and  Knight,  D.,  "Evaluation  of  the  Baldwin-Lomax 

Turbulence  Model  for  Two-Dimensional  Shock  Wave  Boundary  Layer 
Interactions",  AIAA  Paper  No.  83-1697,  AIAA  16th  Fluid  and 
Plasma  Dynamics  Conference,  Danvsers,  Mass.,  July  12-14,  1983. 
Published  in  the  AIAA  J. .  Vol.  22,  July  1984,  pp.  921-928.  [*] 

3.  1  October  1983  -  30  September  1984 

a.  Knight,  D.,  "Numerical  Simulation  of  Three-Dimensional  Shock- 

Turbulent  Boundary  Layer  Interaction  Generated  by  a  Sharp  Fin", 
AIAA  Paper  No.  84-1559,  AIAA  17th  Fluid  Dynamics,  Plasmady- 
namics  and  Lasers  Conference,  June  25-27,  1984.  Published  in 
the  AIAA  vL,  Vol.  23,  December  1985,  pp.  1885-1891.  [*] 

b.  York,  B.,  and  Knight,  D.,  "Calculation  of  TWo-Dimensional  Tur¬ 
bulent  Boundary  Layers  Using  the  Baldwin-Lomax  Model",  AIAA 
23rd  Aerospace  Sciences  Meeting,  Jan  14-17,  1984.  Published  in 
the  AIAA  J..  Vol.  23,  Dec  1985,  pp.  1849-1850. 

4.  1  October  1984  -  30  September  1985 

• 

a.  Knight,  D.,  "Modelling  of  Three  Dimensional  Shock  Wave  Turbu¬ 
lent  Boundary  Layer  Interactions",  in  Macroscopic  Model  ling  of 
Turbulent  Flows.  Lecture  Notes  in  Physics,  Vol  230,  Springer- 
Verlag,  NY,  1985,  pp.  177-201. 


* 
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b.  Knight,  D.,  Horstman,  C.,  Shapey,  B. ,  and  Bogdonoff,  S.,  "The 
Flowfield  Structure  of  the  3-D  Shock  Wave  -  Boundary  Layer 
Interaction  Generated  by  a  20  deg  Sharp  Fin  at  Mach  3",  AIAA 
Paper  No.  86-0343,  AIAA  24th  Aerospace  Sciences  Meeting,  Jan¬ 
uary  6-9,  1986.  Submitted  for  publication  in  the  AIAA  J^_ 

c.  Ong,  C.,  and  Knight,  D.,  "A  Comparative  Study  of 
the  Hybrid  MacCormack  and  Implicit  Beam-Warming  Algorithms  for 
a  Two-Dimensional  Supersonic  Compression  Corner",  AIAA  Paper 
No.  86-0204,  AIAA  24th  Aerospace  Sciences  Meeting,  January  6-9, 
1986.  Submitted  for  publication  in  the  AIAA  J^ 


B.  Interactions  with  Research  Group  at  Princeton  Gas 


cs 


Laboratory  -  1  October  1984  to  1  November  1985 


1.  11  October  1984  :  Meeting  with  Princeton  Gas  Dynamics  Lab 

Research  Group 

Topics  :  1)  Discussion  of  computed  results  for  3-D  sharp  fin  (a  =  20 

deg  and  8  =  0.5  inch). 

2)  Discussion  of  planned  boundary  layer  profile  measurements 
for  the  3-D  sharp  fin  (a  =  20  deg,  S  =  0.5  inch)  at 
Princeton  Gas  Dynamics  Laboratory. 

3)  Discussion  of  C.  Horstman's  calculations  for  the  3-D  swept 
compression  corner. 

4)  Discussion  of  future  collaborative  computational  and 
experimental  research  on  3D  turbulent  interactions 

2.  3  December  1984  :  Conversation  with  S.  Bogdonoff 

Topics  :  1)  Discussion  of  joint  paper  with  Princeton,  Rutgers  and  NASA 

Ames  on  flowfield  structure  of  3-D  sharp  fin 

3.  3  December  1984  :  Conversation  with  S.  Goodwin 

Topics  :  1)  Discussion  of  experimental  surface  pressure  data  for  3-D 

sharp  fin  (a  =  20  deg,  8  =  0.5  inch). 

4.  21  December  1984  :  Conversation  with  S.  Bogdonoff 

Topics  :  1)  Discussion  of  boundary  layer  profile  measurements  for  3-D 

sharp  fin  (a  =  20  deg,  8  =  0.5  inch). 

2)  Discussion  of  future  experimental  measurements  for  3-D 
swept  compression  corner. 


» 


20  February  1985  :  Conversation  with  S.  Bogdonoff 

Topics  :  1)  Discussion  of  initial  experimental  boundary  layer 

measurements  for  3-D  sharp  fin  (a  =  20  deg,  8  =  0.5  inch). 

2)  Discussion  of  recent  measurements  of  3-D  sharp  fin  (<*  = 
17.25  deg) 

22  April  1985  :  Conversation  with  S.  Bogdonoff 

Topics  :  1)  Discussion  of  initial  experimental  boundary  layer 

measurements  for  3-D  sharp  fin  (a  =  20  deg,  8  =  0.5  inch). 

2)  Discussion  of  planned  measurements  for  3-D  swept 
compression  corner. 

21  May  1985  :  Conversation  with  S.  Bogdonoff 

Topics  :  1)  Discussion  of  planned  second  experimental  data  set  for  3-D 

sharp  fin  (a  =  20  deg,  8  =  0.5  inch). 

2)  Discussion  of  planned  computation  of  3-D  swept  compression 
corner  for  (oc,X)  =  (24,60)  deg. 

3)  Discussion  of  experimental  investigation  of  17.25  deg 
sharp  fin,  25  deg  semi-cone  and  (oc.x)  =  (30,60)  swept 
compression  corner. 

1  August  1985  :  Conversation  with  S.  Bogdonoff 

Topics  :  1)  Discussion  of  structure  of  3-D  sharp  fin  flowfield. 

22  August  1985  :  Conversation  with  S.  Bogdonoff 

Topics  :  1)  Discussion  of  flowfield  structure  for  3-D  swept 


compression  corner. 


2)  Discussion  of  planned  experimental  investigation  of 
symmetric  3-D  sharp  fin. 

3)  Discussion  of  future  computations  for  3-D  swept 
compression  corner. 

10.  23  September  1985  :  Conversation  with  S.  Bogdonoff 

Topics  :  1)  Discussion  of  3-D  turbulent  interaction 

flowfield  structures  for  sharp  fin  and  swept  compression 
corner . 

11.  9  Oct  1985  :  Meeting  with  Princeton  Gas  Dynamics  Lab 

Research  Group 

Topics  :  1)  Discussion  of  experimental  surface  visualization  (kerosene 

lampblack)  for  3-D  sharp  fin 

2)  Discussion  of  experimental  schlieren  photographs  for  3-D 
sharp  fin 

12.  21  Oct  1985  :  Meeting  with  Princeton  Gas  Dynamics  Lab 

Research  Group,  C.  Horstman  (NASA  Ames) 

Topics  :  1)  Discussion  of  computed  particle  pathlines  for  3-D  sharp 

fin  (a  =  20  deg,  8  =  0.5  inch). 

2)  View  videotape  of  particle  pathlines 

2)  Development  of  model  for  mean  flowfield  structure  for  3-D 
sharp  fin 
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Interactions  with  C.  C.  Horstman  (NftSA  Ames  Research  Center] 

1  October  1984  -  1  November  1985 

16  Oct  1984  :  Conversation  with  Mike  Horstman 

Topics  :  1)  Current  computations  of  3-D  sharp  fin  for  fin  angles  of  10 

to  20  deg 

2)  Discussion  of  secondary  separation  line 
structure  for  3-D  sharp  fin 

28-30  November  1984  :  Conversations  with  Mike  Horstman 

Topics  :  1)  Discussion  of  status  of  computations  for  3-D  sharp  fin  by 

Horstman  and  Knight 

2)  Discussion  of  future  computations 

4  February  1985  :  Conversation  with  Mike  Horstman 

Topics  :  1)  Discussion  of  comparison  of  computed  results  by  Knight  and 

Horstman  for  3-D  sharp  fin  (a  =  20  deg,  8  =  0.5  inch) 

14  February  1985  :  Conversation  with  Mike  Horstman 

Topics  :  1)  Further  discussion  of  computed  results  for  3D  sharp  fin. 

2)  Requested  additional  computed  results  from  Horstman  fgr 
comparison 

3)  Discussion  of  future  work 

28  February  1985  :  Conversation  with  Mike  Horstman 

Topics  :  1)  New  wall  function  turbulence  model  employed  by  Horstman 


25  March  1985  :  Conversation  with  Mike  Horstman 

Topics  :  1)  Effect  of  new  wall  function  model  on  predicted  results  for 

swept  compression  corner 

16  April  1985  :  Visit  by  Mike  Horstman  to  Rutgers  University 

Topics  :  1)  Flowfield  structure  for  3-D  sharp  fin  interaction 

2)  Contents  of  joint  paper  with  Knight,  Horstman,  Bogdonoff 
and  Shapey  for  AIAA  24th  Aerospace  Sciences  Meeting  in  Jan 
1986 

3)  Discussion  of  recent  computated  results  for  3-D  swept 
compression  corner 

4)  Discussion  of  2-D  compression  ramp,  including  calculations 
by  C.  Ong  (Rutgers),  and  Dennis  Johnson  (NASA  Ames). 

5)  Discussion  of  experimental  fluid  dynamics  program  at  NASA 
Ames 

25  April  1985  :  Conversation  with  Mike  Horstman 

Topics  :  1)  Contents  of  joint  paper  with  Knight,  Horstman,  Bogdonoff 

and  Shapey  for  AIAA  24th  Aerospace  Sciences  Meeting  in  Jan 
1986 

2)  Discussion  of  computational  results  for  3-D  swept 
compression  corner  (oc.X)  =  (24,60)  deg 

10  and  19  June  1985  :  Conversation  with  Mike  Horstman 

Topics  :  1)  Discussion  of  new  3-D  finite  volume  code  being  developed 

by  Mike  Horstman 


2)  Discussion  of  parameters  for  3-D  sharp  fin  calculations 


10.  16  July  1985  :  Conversation  with  Mike  Horstman 

Topics  :  1)  Boundary  conditions  employed  by  Horstman  for  3-D  swept 

compression  corner  calculations 

11.  16  August  1985  .-  Cornmunication  with  Mike  Horstman 

Topics  :  1)  Sent  Mike  Horstman  the  results  of  computations  by  Knight 

for  3-D  swept  compression  corner  (ot.x)  =  (24,60)  deg 

12.  22  August  1985  :  Cornmunication  with  Mike  Horstman 

Topics  :  1)  Sent  Mike  Horstman  additional  results  of  computations  for 

3-D  swept  compression,  with  emphasis  on  effect  of  upstream 
boundary  conditions 

13.  6  September  1985  .-  Conversation  with  Mike  Horstman 

Topics  :  1)  Discussion  of  comparison  of  computed  (Knight  and  Horstman) 

and  experimental  (Princeton  Gas  Dynamics  Lab)  data  for  3-D 
sharp  fin  (a  =  20  deg,  8  =  0.5  inch) 

14.  16  Oct  1985  .-  Communication  with  Mike  Horstman 

Topics  .•  1)  Sent  comparison  of  experimental  and  computed  profiles  of 

pitot  pressure  and  yaw  angle  for  3-D  sharp  fin  (<x  =  20 
deg,  8  =  0.5  inch).  Copies  also  sent  to  Shapey  and 
Bogdonoff  (Princeton). 

15.  22  Oct  1985  :  Visit  by  Mike  Horstman  to  Rutgers  University 

Topics  :  1)  Dis<russion  of  particle  pathlines  for  3-D  sharp  fin 

2)  Discussion  of  flowfield  structure  for  3-D  sharp  fin 


Discussion  of  videotape  of  particle  pathlines  to  be 
presented  at  AIM  Aerospace  Sciences  Meeting  in  Jan  1986 

Discussion  of  contents  of  paper  by  Knight,  Horstman, 
Shapey  and  Bogdonoff  to  be  presented  at  AIAA  24th 
Aerospace  Sciences  Meeting  in  Jan  1986 

Discussion  of  new  concepts  in  3-D  flowfield  graphics 

Discussion  of  flowfield  structure  of  3-D  swept  compression 
corner 
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D.  Spoken  Papers  Presented  at  Technical  Meetings 
1  October  1984  to  1  November  1985 

1.  Knight,  D.,  "Calculation  of  Three-Dimensional  Shock  Wave-Turbulent 
Boundary  Layer  Interaction  Generated  by  a  Sharp  Fin",  Thirty-Seventh 
Annual  Meeting,  Division  of  Fluid  Dynamics,  American  Physical  Society, 
18-20  November  1984,  Bulletin  of  American  Physical  Society,  Vol.  29,  No. 
9,  Nov  1984,  p.  1569. 

2.  Knight,  D.,  "Theoretical  Investigation  of  Shock  Wave  -  Turbulent 
Boundary  Layer  Interaction  :  Problems  in  Reconciling  Computation  and 
Experiment",  Workshop  on  Structure  of  High  Speed  Turbulent  Boundary 
Layers,  Princeton  University,  29-30  July  1985. 

E.  Seminars  -  1  October  1984  to  1  November  1985 

1.  Knight,  D.,  "Theoretical  Investigation  of  3-D  Shock  Wave-Turbulent 
Boundary  Layer  Interaction  Generated  by  a  Sharp  Fin",  Air  Force  Wright 
Aeronautical  Laboratory,  Wright-Patterson  AFB,  OH,  28  October  1985. 
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Fig.  3  Comparison  of  Computed  and  Experimental  Surface  Pressure 
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Abstract 

A  comparative  study  Is  made  between  the 
MacCormack  explicit-implicit  predictor-corrector 
and  the  Beam-Warming  fully  implicit  algorithms  for 
solving  compressible  viscous  flow.  The  mass- 
averaged  two-dimensional  compressible  Navler-Stokes 
equations  In  strong  conservation  law  form  and 
general  curvilinear  coordinates  are  solved 
numerically  by  marching  forth  In  time  on  a 
body-fitted  curvilinear  grid  for  a  shock  wave- 
turbulent  boundary  layer  Interaction  over  a 
two-dimensional  compression  corner.  Along  the 
surface  of  the  corner  the  boundary  condition  for 
the  implicit  part  of  the  hybrid  MacCormack 
algorithm  Is  formulated  using  an  approximation 
Involving  a  lag  in  time  of  one-half  time  step. 
Turbulence  Is  simulated  by  means  of  the 
Baldwln-Lomax  algebraic  turbulent  eddy  viscosity 
model.  Computations  are  performed  for  a  Mach 
number  of  1.96  with  a  Reynolds  number  Re4  (based 
on  the  Incoming  boundary  layer  thickness  <•)  of 
0.25  x  10°,  and  for  a  Mach  number  of  2.83  with  a 
Reynolds  number  of  1.8  x  10®.  The  primary 
objectives  of  the  study  are,  1)  to  determine  the 
extent  to  which  the  steady  state  solution  obtained 
by  the  hybrid  MacCormack  algorithm  Is  dependent 
upon  the  size  of  the  time  step  employed  in  marching 
the  calculation  toward  the  steady  state  solution, 
2)  to  compare  the  two  algorithms  regarding  accuracy 
and  efficiency,  and  3)  to  further  examine  the 
efficacy  of  the  Baldwln-Lomax  algebraic  turbulent 
eddy  viscosity  model  through  comparison  with  recent 
experimental  measurements  of  the  Reynolds  shear 
stress. 


I.  Introduction 

Until  the  mid  1970' s  most  numerical  schemes 
used  to  solve  the  Navler-Stokes  equations  were 
explicit  algorithms  such  as  MacCormack's  explicit 
method*.  The  restriction  placed  on  the  size  of  the 
maximum  allowable  time  step  by  the  CFL  (Courant, 
Lewy,  and  Friedrichs)  condition  severely  limited 
the  usefulness  of  explicit  procedures  In  the 
solution  of  high  Reynolds  number,  viscous, 
compressible  flow.  8y  the  mid  1970' s  this  CFL 
restriction  was  removed  by  the  fully  explicit 
numerical  schemes  of  Briley  and  McOonald*,  and  Beam 
and  Warming^  which  became  very  widely  used.  In 
particular,  the  8eam-Wanm1ng®  method  has  been  used 
to  compute  flows  around  airfoils*,  nosetlps®, 
boattall  afterbodies*,  cascades7,  and  Inside  2-0 
inlets*  and  nonaxlsymmetrlc  nozzles*.  However, 
these  Implicit  methods  require  linearization  of 
terms  In  the  governing  equations  In  order  to  form 

'Graduate  Student,  Oept  of  Mechanical  and  Aero¬ 
space  Engineering;  Member,  AIAA.  Presently  at 
Continuum  Dynamics  Inc.,  Princeton,  NJ. 
**Professor,  Oept.  of  Mechanical  and  Aero¬ 
space  Engineering;  Associate  Fellow,  AIAA. 


desirable  block  matrix  structures.  The  resulting 
block  pentadlagonal  or  septadlagonal  matrices  need 
to  be  approximately  factored  Into  block  tridlagonal 
matrices  before  an  efficient  matrix  Inversion 
procedure  can  be  applied.  Both  these 
approximations  result  In  limitations  to  the  size  of 
the  maximum  allowable  time  step*®. 

In  1981  MacCormack**  presented  an 
explicit-implicit  predictor-corrector  method  whicn 
Involved  the  simple  inversion  of  block  bidiagonal 
matrices  In  an  effort  to  further  reduce  demand  on 
computer  time.  Since  Its  introduction  the  hybr‘d 
scheme  has  been  used  by  numerous  investigators  to 
compute  a  large  variety  of  compressible  viscous 
flows.  In  calculating  a  separated  shock  wave- 
laminar  boundary  layer  interaction  over  a  flat 
plate  MacCormack  was  able  to  reduce  the  required 
computer  time  by  a  factor  of  17.5  relative  to  his 
fully  explicit  method  while  maintaining  comparable 
accuracy.  Shang  and  MacCormack*7  evaluated  the  new 
method  against  Its  fully  explicit  predecessor  for  a 
Mach  8  flow  over  an  axlsymmetric  biconic  body  and 
achieved  computer  time  reduction  by  a  factor  of  13. 

Accuracy  and  efficiency  are  but  two  features 
desirable  In  a  numerical  method.  Independence  of 
the  steady  state  solution  from  the  size  of  the 
chosen  time  step  Is  another.  Before  asking  whether 
the  Steady  state  solution  of  a  given  numerical 
algorithm  Is  accurate  it  seems  logical  to  first 
Inquire  If  it  yields  the  same  steady  state  solution 
Irregardless  of  the  time  step  size  selected. 
MacCormack  advocated****®  a  successive  reduction  of 
the  time  step  size  near  the  end  of  the  calculation 
of  high  Reynolds  number  flows  In  order  to  avoid  any 
possible  dependence  of  the  steady  state  solution  on 
the  time  step  size.  This  concern  about  time  step 
dependence  was  underscored  by  Kumar**  who  reported 
finding  considerable  time  step  dependence  of  his 
steady  state  solution  near  the  immediate 
neighborhood  of  a  separation  region  induced  by 
shock  Impingement  upon  the  turbulent  boundary  layer 
In  a  duct  whose  inflow  Mach  number  was  5.  Howevgr, 
no  significant  CFL  dependence  was  found  by  Gupta. 
Gnoffo,  and  MacCormack**  who  computed  an 
unseparated  laminar  flow  over  an  ax i symmetric  at 
Mach  number  44  with  a  detached  bow  shock. 

Other  Investigators  who  used  MacCormack's 
expl Id t-lmol Id t  scheme  Included  Kordulla  and 
MacCormack**;  White  and  Anderson*7;  Hung  and 
Kordulla*®:  and  Imlay,  Kao,  McMaster,  and 
MacCormack*®.  Despite  the  indication  of  possible 
significant  CFL  dependence  of  the  steady  state 
solution  by  MacCormack  and  Kumar,  few  Investigators 
have  addressed  this  issue  and  some  have  stated 
without  further  substantiation  that  "the  implicit 
MacCormack  method  Is  unreliable  such  that  the 
steady  state  solutions  depend  on  time 
increments”7®.  Nevertheless,  it  is  obviously 
important  to  ascertain  the  severity  of  any  such 
dependence,  particularly  in  the  presence  of 


separated  regions.  Hence,  the  first  objective  of 
the  present  study  fs  to  determine  whether  the 
steady  state  solution  obtained  by  the  MacCormack 
explicit-implicit  method  depends  on  the  time  step 


Since  the  8eam-Warm1ng  method  has  become  very 
popular  and  the  explicit-implicit  hybrid  method 
appears  to  be  accurate  and  efficient,  the  question 
regarding  their  relative  accuracy  and  efficiency  Is 
clearly  a  relevant  and  Important  one.  Lawrence, 
Tannehlll,  and  Chausee?!  discussed  this  Issue  In 
the  case  of  the  Mach  2  laminar  flow  over  a  flat 
plate.  They  used  both  algorithms  to  solve  the 
parabolized  Navler-Stokes  equations  by  marching  In 
space  rather  than  In  time  for  the  Beam-Harming  and 
MacCormack  hybrid  methods.  They  observed  that  the 
computer  time  requirements  for  the  Beam-Harming  and 
MacCormack  hybrid  were  approximately  equal.  Iyer 
and  von  Lavante“  attempted  a  comparison  for  a 
viscous  transonic  flow  in  turbomachinery  cascasdes. 
However,  a  more  extensive  study  for  a  time-marching 
solution  of  the  full  two-dimensional  Navler-Stokes 
equations  over  a  realistic  grid  appears  necessary. 
Hence,  the  second  objective  of  the  present  study  Is 
an  accuracy  and  efficiency  comparison  between  the 
MacCormack  explicit-implicit  and  the  Beam-Harming 
fully  Implicit  algorithms. 

The  configuration  selected  for  this 
comparative  study  Is  supersonic  turbulent  flow  past 
a  2-0  compression  corner  (Fig.  1).  Related  flows 
are  common  In  many  engineering  applications 
Including  turbomachinery  and  high  speed  aircraft^. 
In  view  of  the  practical  Importance  of  shock 
wave-turbulent  boundary  layer  Interactions,  It  Is 
not  surprising  therefore  that  this  configuration 
has  been  extensively  Investigated  both 
experimentally  and  theoretically.  Several  reviews 
have  been  published.  Including  Green^*,  Korkegl’3, 
and  Hankey  and  Holden25.  jn  particular,  extensive 
experimental  measurements  have  been  performed  by 
Settles’®  for  the  specific  case  of  a  Mach  2.83 
turbulent  flow  past  a  2-0  compression  corner  at  a 
series  of  corner  angles  o  *  8  to  24  deg,  and  for  a 
range  of  Reynolds  nunbers  Res^  =  0.76  x  10®  to  7.7 
x  10®.  Settles'  measurements  Include  surface 
properties  (pressure,  skin  friction  and  surface  oil 
flow  visualization),  and  boundary  layer  profiles  of 
velocity,  static  pressure  and  Mach  number. 
V1sbal”»28  computed  the  entire  set  of  flows 
corresponding  to  the  experimental  configurations  of 
Settles,  utilizing  the  algebraic  turbulent  eddy 
viscosity  model  of  Baldwin  and  Lomax.  Vlsbal 
observed  1)  the  8aldw1n-Lomax  outer  function  was 
unsuitable  for  determination  of  the  length  scale  of 
the  turbulence  In  the  separation  region,  2) 

1  ncorporatlon  of  the  relaxation  model  of  Shang  and 
Hankey29  improved  the  prediction  of  the  extent  of 
upstream  propagation  of  the  disturbance  associated 
with  the  corner,  and  3)  the  computed  boundary  layer 
recovery,  downstream  of  reattachment,  was 
significantly  less  than  observed  experimentally. 

Recently,  measurements  of  turbulent  Reynolds 
stresses  have  been  performed  for  the  same 
configuration  by  Hayakawa,  Muck,  Smlts  and 
Sp1na30,31,32.  These  measurements,  therefore, 
provide  the  opportunity  for  direct  examination  of 
the  efficacy  of  the  8aldw1n-Lomax  turbulence  model. 
Hence,  the  third  objective  of  the  present  study  Is 
to  directly  compare  the  computed  and  measured 
Reynolds  shear  stress  for  the  2-0  compression 


corner  at  Mach  2.B3  for  several  corner  angles,  and 
to  attempt  to  elucidate  the  deficiencies  In  the 
Saldwln-Lomax  model.  In  addition,  computations 
have  been  performed  for  2-0  compression  corner  at 
Mach  1.96  and  Re$  *  0.25  x  10®  for  o  *  16  deg,  ana 
the  results  compared  with  the  experimental  data  of 
0oll1ng33  for  surface  pressure. 


II.  Method  of  Solution 


Governing  Equations 


The  flow  was  assumed  to  be  described  by  the 
two-dimensional,  mass-averaged  compressible 
Navler-Stokes  equations.  The  equations  were 
written  In  strong  conservation  law  form  and  general 
curvilinear  coord1nates3*.35.  the  fluid  was 
assumed  to  be  a  perfect  gas  with  a  molecular 
dynamic  viscosity  given  by  Sutherland's  law  and  a 
constant  molecular  Pranatl  number  of  0.72. 


Turbulence  was  simulated  by  the  Baldwln- 
Lomax36  algebraic  eddy  viscosity  model  with  r  = 
0.40,  A*  *  26,  C|(ie|3  =  0.3  and  k  =  0.0168.  The 
value  of  CCp  was  obtained  from  the  investigation  of 
fork  and  Knight^'  and  Vlsbal^  wh1ch  detailed  the 
dependence  of  Ccp  on  Mach  number.  In  particular, 
CCn  varies  from  1.8  at  Mach  2  to  2.1  at  Mach  3.  The 
relaxation  model  of  Shang  and  Hankey^®  was  employed 
In  all  cases  except  the  8  deg  compression  ramp  at 
Mach  2.83,  with  the  relaxation  length  scale  set 
equal  to  the  upstream  boundary  layer  thickness.  The 
turbulent  Prandtl  number  Is  0.9. 


putatlonal  Domain  and  Boundary  Conditions 


Nearly  orthogonal,  body-fitted  grids  were 
generated  using  the  method  of  Vlsbal  and  Knight®. 
Typically,  30  grid  points  were  contained  within  the 
Incoming  boundary  layer,  and  the  distance  between 
the  wall  and  the  nearest  line  of  grid  points  had  a 
T*  value  smaller  than  2.5,  where  Y+  Is  the  distance 
normal  to  the  wall  nond 1 mens  1 ona 1 1  zed  by  the  local 
friction  velocity  and  the  wall  value  of  the 
kinematic  viscosity.  For  each  flow  configuration, 
the  same  grid  was  employed  for  the  calculations 
using  the  Beam-Harming  and  MacCormack  hybrid 
algorithm. 


The  inflow  boundary  was  positioned  in  the 
undisturbed  turbulent  flat  plate  boundary  layer 
where  the  computed  momentum  thickness  matched  the 
experimental  value.  In  both  the  Beam-Warming  and 
hybrid  MacCormack  computations,  the  flow  variables 
on  the  Inflow  boundary  were  held  fixed  at  the  given 
values.  The  Inflow  boundary  condition  for  the 
Implicit  step  of  the  Hybrid  MacCormack  was 
prescribed  by  setting  the  temporal  change  in  the 
solution^*,  4U,  to  zero.  The  outflow  boundary 
was  located  far  enough  downstream  of  the  corner  to 
be  In  a  region  of  small  streamwlse  flow  gradients. 
Along  the  outflow  boundary  the  extrapolation 
condition  was  assumed.  For  this  boundary,  both  the 
Beam-Warming  and  the  explicit  step  of  the  hybrid 
MacCormack  represented  JU/3C  *  0  by  a  first  order 
accurate  dl  fferencl ng,  where  C  Is  the  transformed 
coordinate  In  the  general  streamwlse  direction  and 
U  Is  the  vector  of  dependent  variables 3*11.  The 
Implicit  part  of  MacCormack  s  method  set  the 
temporal  change  ail  at  the  downstream  boundary 
equal  to  Its  value  at  the  adjacent  constant-C  line. 


At  the  lower  boundary,  the  velocity  and  normal 
derivative  of  the  static  pressure  were  set  to  zero. 
For  the  Mach  1.96  compression  corner  flow, 
adiabatic  boundary  conditions  were  used  for  both 
schemes,  while  a  constant  (near  adiabatic) 
temperature  was  specified  for  the  Mach  2.83  corner. 
The  lower  boundary  condition  for  the  Implicit 
portion  of  MacCormack's  hybrid  was  formulated  by 
allowing  a  one-half  time  step  lag  in  the  value  of 
the  temporal  change  atJ.  The  upper  boundary  was 
placed  sufficiently  far  from  the  lower  so  that 
freestream  conditions  prevailed  all  along  its 
length.  For  the  Beam-Warming  scheme,  a  non- 
reflection  condition*?  was  applied  there.  For  the 
hybrid  MacCormack,  however,  the  flow  variables 
along  the  upper  boundai^  were  simply  held  at  the 
freestream  value  and  4U  was  set  to  zero.  The 
shock  emerges  through  the  outflow  boundary. 

Numerical  Procedure 

The  fully  implicit  scheme  studied  was  the 
approximate  factori zati on  algorithm  of  Beam- 
Warming?,  formulated  using  Euler  implicit 
time-differencing  and  second-order  accurate, 
centered  differencing  for  the  spatial  derivatives. 
Fourth  order  explicit  damping  terms  were  Included 
in  the  manner  shown  by  Thomas". 

The  expl ici t-lmpl 1  cl t  algorithm  was  that  of 
MacCormackfl ,  extended  to  2-0  general  coordinates 
by  von  Lavante  and  Thompkins^.  The  algorithm  is 
second  order  accurate  in  space  and  time.  While 
marching  in  time  the  order  of  finite  differencing 
in  the  explicit  predictor  and  corrector  steps  was 
cycled  from  one  step  to  the  next  while  that  in  the 
implicit  steps  was  kept  as  forward  differencing  in 
the  predictor  and  backward  In  the  corrector.  At 
all  times  opposite  orders  of  differencing  were 
employed  In  the  predictor  and  corrector.  The  usual 
fourth  order  damping  expressed  in  terms  of  the 
pressure  is  used  for  the  explicit  part.  Implicit 
damping  was  also  incorporated,  in  the  manner 
suggested  by  MacCormackH. 

Both  the  explicl t-impl ici t^O  and  fully 
Implicit*?  computer  codes  were  carefully  validated 
with  excellent  accuracy  for  a  variety  of  flows, 
including  laminar  and  turbulent  boundary  layers, 
and  shock-laminar  boundary  layer  interaction. 


III.  Results  and  Discussion 
Courant  Number  Dependence  of  the  Hybrid  MacCormack 

As  indicated  above,  the  first  objective  of  the 
research  is  to  examine  the  possible  Courant  number 
dependence  of  the  steady  state  solution  computed  by 
the  hybrid  MacCormack  algorithm.  During  this 
examination  it  will  be  convenient  to  also  examine 
the  accuracy  of  the  solution  in  comparison  with 
Beam-Warming  results  and  the  experimental 
measurements  of  0o111ng33  and  Settles  et  al.**. 

Mach  1.96  flow  over  16  deg  compression  corner 

The  computed  and  measured  surface 
pressure  distributions  for  the  Mach  1.96 
compression  corner  are  displayed  in  Fig.  2.  The 
Reynolds  number  Res  *  0.25  x  IQ**  and  a  *  16  deg. 
In  this  and  all  subsequent  figures,  X  denotes  the 
distance  from  the  corner  measured  along  the 


surface.  Calculated  profiles  shown  are  for  the 
MacCormack  hybrid  method  at  Courant  numoers  of  0.9 
(fully  explicit)  and  45  (hybria),  and  the 
Beam-Warming  algorithm,  where  the  Courant  numoer  is 
defined  by  Shang^V  The  experimental  data  of 
0olling33  4re  also  shown.  The  results  clearly 
indicate  that  the  steady  state  solution  for  the 
surface  pressure  using  the  MacCormack  explicit- 
implicit  algorithm  Is  insensitive  to  the  Courant 
number  and  very  close  to  the  Beam-Warming  results. 
The  computed  upstream  propagation  of  the  surface 
pressure,  measured  from  the  corner  (X  *  0),  is 
approximately  301  below  the  experimental  value. 
Since  the  extent  of  the  upstream  propagation  is 
directly  related  to  the  magnitude  of  the  length 
scale  employed  in  the  relaxation  odel,  the  length 
scale  4.  Is  too  small  for  this  ase.  The  present 
results  at  Mach  2  for  Re5  *  0.25  x  10&,  together 
with  previous  results  ofl/isbal*®  at  Mach  3  for 
ReJa>  *  0.76  x  106  to  7.7  x  10^,  imply  that  the 
relaxation  length  is  a  moderate  function  of  Re4^ 
(i.e.,  the  relaxation  length*?  increases  witn 
decreasing  Re4  ).  This  observation  is  consistent 
with  the  resuTts  of  Shang  and  Hankey*^,  who 
employed  a  relaxation  length  of  104m  for  their 
studies  of  the  2-0  compression  corner  at  Mach  3  for 
Res-  ■  0.14  x  106. 

The  calculated  skin  friction  coefficient  Cf 
distribution  for  the  same  flow  is  shown  in  Fig.  3. 
The  computed  results  using  the  MacCormack  hybrid 
algorithm  are  again  observed  to  have  no  marked 
dependence  on  the  Courant  number  despite  the  fact 
that  the  two  Courant  numbers  differ  by  a  factor  of 
50.  The  computed  results  are  in  good  agreement 
with  the  computation  using  Beam-Warminq  s  fully 
Implicit  algorithm.  The  results  using  MacCormack's 
method  manifest  a  small  streamwise  oscillation  in 
Cf  downstream  of  reattachment.  The  cause  of  this 
oscillation  is  currently  under  Investigation. 

In  Fig.  4,  the  computed  velocity  parallel  to 
the  wall,  normalized  by  the  upstream  freestream 
velocity  U»,  is  displayed  at  X/4„  =  0.16  (down¬ 
stream  of  the  corner)  for  the  Beam-Warming  and 
MacCormack  hybrid  algorithms.  In  the  latter  case, 
results  are  displayed  for  Courant  numbers  of  0.9 
and  45.  In  this  and  all  subsequent  figures,  the 
distance  normal  to  the  surface  is  denoted  by  y .  The 
results  indicate  that  the  computed  solution  using 
the  MacCormack  hybrid  algorithm  is  insensitive  to 
the  Courant  number,  and  in  close  agreement  with  the 
results  obtained  using  the  Beam-Warming  algorithm. 
Similar  conclusions  were  obtained  by  examination *of 
static  temperature,  Baldwin-iomax  outer  function, 
static  pressure  and  eddy  viscosl  ty<0. 

Mach  2.83  flow  over  16  deg  compression  corner 

The  surface  pressure  distributions  for  a  Mach 
2.83  flow  at  Re4  *  1.6  x  10*>  and  a  =  16  deg  are 
shown  in  Fig.  5  f^r  computations  using  the  methods 
of  MacCormack  and  Beam-Warming.  The  experimental 
measurements  of  Settles  et  a  1 . 26  are  also 
displayed.  A  high  Courant  number  of  85  and  a  low 
of  30  used  by  the  hybrid  algorithm  give  no 
discernable  difference  in  the  computed  surface 
pressure.  There  is  also  good  agreement  with  the 
pressure  calculated  using  the  Beam-Warming  scheme 
as  well  as  with  experiment. 

In  Fig.  6  the  corresponding  skin  friction 
distributions  are  exhibited.  As  before,  the  skin 
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Table  1.  Comouter  Time  Comparison,  Mach 
1.96  Flow  Over  16  deg  Ramp 


friction  computed  by  the  hybrid  method  is 
Insensitive  to  the  time  step  size.  There  is  a 
slight  tendency  for  the  higher  CFl  case  to  predict 
a  marginally  lower  skin  friction  further  downstream 
from  reattachment .  The  hybrid  method  predicts  a 
modestly  higher  skin  friction  than  the  8eam-Warmlng 
method  further  downstream  from  reattachment  but 
agrees  closely  with  the  latter  practically 
everywhere  else. 

In  Fig.  7,  the  computed  and  experimental 
horizontal  velocity  profiles  along  a  vertical  line 
are  displayed  at  the  corner.  Again,  the  results 
clearly  display  no  Courant  number  dependence  for 
the  MacCormack  hybrid  algorithm,  and  are  in  close 
agreement  with  the  profile  calculated  using  the 
Beam-Warming  algorithm.  The  computed  profiles  are 
in  good  agreement  with  experiment^.®  except  In 
the  Itimedlate  vicinity  of  the  surface. 

Mach  2.83  flow  over  20  deg  compression  corner 

The  surface  pressure  for  the  Mach  2.B3  flow  at 
Re4^  *  1.6  x  10^  and  a  ■  20  deg  Is  detailed  In  Fig. 
8.  "This  ramp  angle  Is  the  largest  studied  In  this 
Investigation.  The  Courant  numbers  used  for  the 
McCormack  scheme  are  30  and  70.  The  Beam-Warming 
and  experimental  results®  are  also  plotted.  As  In 
the  previous  cases,  the  computed  surface  pressure 
exhibits  no  Courant  number  dependence  for  the 
MacCormack  algorithm.  Also,  the  results  obtained 
from  the  Beam-Warming  and  MacCormack  hybrid 
algorithms  are  In  close  agreement. 

In  Fig.  9,  the  skin  friction  for  the  same  flow 
Is  shown  for  both  the  MacCormack  method,  the 
Beam-Warming  method,  and  the  experiment.  The 
resu'ts  indicate  that,  even  In  the  presence  of  such 
a  strong  adverse  pressure  gradient  and  large 
separation  region,  the  computed  skin  friction  Is 
Insensitive  to  time  step  size.  The  hybrid  method 
predicts  a  slightly  higher  skin  friction  than  the 
Beam-Warming  scheme  downstream  of  reattachment. 

Efficiency  of  the  Hybrid  MacCormack  Algorithm 

In  addition  to  Courant  number  dependence  and 
accuracy,  the  efficiency  of  the  hybrid  MacCormack 
algorithm  was  also  examined  for  the  Mach  1.96  flow 
over  a  16  deg  ramp.  Computations  using  the  hybrid 
scheme  were  made  at  Courant  nunders  of  0.9  (fully 
explicit),  5.,  10.,  21.,  40.,  and  45.  One 
computation  using  the  Beam-Warming  scheme  was 
performed  which  employed  a  maxlmun  Courant  number 
of  33.  In  order  to  avoid  nvanerlcal  Instability,  It 
was  necessary  to  start  the  8eam-Warmlng  calculation 
at  a  smaller  time  step  and  progressively  Increase 
it  to  a  maximum  consistent  with  numerical 
stability.  Both  numerical  codes  were  written  In 
FORTRAN,  and  executed  on  an  NAS  AS/9000  mainframe 
computer.  A  uniform  set  of  convergence  criteria 
was  employed  for  all  calculations.  It  was  observed 
that  convergence  to  steady  state  required 
approximately  the  same  physical  time  of  Integration 
in  all  cases. 

The  comouter  time  requirements  for  these 
computations  are  tabulated  In  Table  1.  It  is 
observed  for  this  case  that  the  hybrid  scheme 
requires  one-third  of  the  comouter  time  used  by  the 
8eam-Warm1ng  method,  and  up  to  a  factor  of  36  less 
than  the  ful ly-expl Id t  MacCormack  algorithm. 


Algorithm 

Courant 

Number*! 

Computer  Time 

In  hours 
(NAS  AS/9000) 

8eam-Warm1ng 

33. 

8.7 

45. 

2.9 

MacCormack 

40. 

3.2 

Explicit- 

21. 

6.0 

Impl Icit 

10. 

11.8 

Hybrid 

5. 

12.4 

0.9 

(expl Icit) 

104.1 

Reynolds  Shear  Stress  Comparison 

Reynolds  shear  stress 

In  the  present  section,  the  computed  ana 
experimental  profiles  for  the  Reynolds  shear 
stress,  defined  as  -»  u'v',  are  displayed.  The 
quantities  u'  and  v'  are  the  temporal  fluctuating 
velocity  components  parallel  and  normal  to  the 
wall,  respectively.  The  overbar  represents  the 
time  average.  The  experimental  Reynolds  shear 
stress  Is  obtained  from  the  measurements^.^?  of 
-(pu)'v'  by  employing  the  "Strong  Reynolds  Analogy" 
(l.e.,  pressure  fluctuations  are  small  compared  to 
density  or  temperature  fluctuations)  and  "Very 
Strong  Reynolds  Analogy"  (l.e.,  fluctuations  in 
total  temperature  are  neglected).  The  uncertainties 
(approximately  r30X)  In  the  measurement  of  the 
Reynolds  stress  are  discussed  In  Ref.  30-32.  The 
theoretical  Reynolds  shear  stress,  modeled  using 
the  Baldwln-Lomax  algebraic  turbulent  eddy 
viscosity.  Is  -pu"v",  where  u"  and  v"  denote  the 
mass-averaged  fluctuating  velocity  components 
parallel  and  normal  to  the  wall,  respectively.  With 
the  assumption  of  the  Strong  and  Very  Strong 
Reynolds  Analogies,  the  theoretical  Reynolds  shear 
stress  Is  approximately  equal  to  -p  u'v'.  In  all 
plots,  the  experimental  and  theoretical  Reynolds 
shear  stress  are  normalized  by  0.5  p«U^. 

Mach  2.83  flow  over  8  deg  compression  corner 

The  computed  and  experimental  Reynolds 
shear  stress  profile  at  x  *  -0.63  located 
upstream  of  the  Interaction  region.  Is  displayed  in 
Fig.  11  for  the  Mach  2.83  flow  over  an  8  degree 
compression  corner.  The  computed  results  of 
MacCormack 's  hybrid  method  at  Courant  numbers  of  95 
and  0.9  are  plotted  together  with  the  Beam-Warming 
and  experimental  profiles  of  Muck  et  al.^O-32. 
Agreement  among  all  three  computed  Reynolds  stress 
profiles  is  excellent.  Agreement  between  the 
computed  and  experimental  profiles  Is  good  exceot 
In  the  outer  portion  of  the  boundary  layer  where 
the  predicted  Reynolds  stress  is  low. 

In  Figs.  12  and  13,  computed  and  experimental 
Reynolds  stress  profiles  are  shown  at  x/j„  *  0.78 
and  1.17.  They  indicate  that  the  predicted  peak 
value  of  the  Reynolds  stress  profile  is 
significantly  too  low,  and  that  the  predicted  peak 
Is  located  too  near  the  wall. 


Mach  2.83  flow  over  16  deg  compression  corner 

For  the  Mach  2.83  flow  over  a  16  deg 
compression  corner,  the  calculated  and  experimental 
Reynolds  stress  profiles  at  stations  %/»„,  =  0.49, 
2.0  and  5.4  are  exhibited  In  Figs.  14  to  16.  It  Is 
apparent  that  the  peaks  of  the  computed  and 
experimental  Reynolds  stress  profiles  are 
comparable  in  magnitude.  This  Is  summarized  In 
Fig.  17,  which  displays  the  distribution  of  the 
magnitude  of  the  peak  of  the  Reynolds  stress 
profiles  with  distance  X.  It  Is  evident,  however, 
that  the  computed  peak  Is  located  too  close  to  the 
wall.  The  height  of  the  peak  Vpeai(  as  a  function  of 
distance  X  Is  displayed  In  Fig.  18.  The  distance 
of  the  peak  from  the  wall  Increases  as  the  boundary 
layer  develops  downstream  for  both  computed  and 
experimental  profiles. 

Mach  2.33  flow  over  20  deg  compression  corner 

[n  Figs.  19  to  21,  the  profiles  'of  the 
Reynolds  shear  stress  at  stations  X/4*  =  1.0,  2.3 
and  4.6,  are  shown  for  the  Mach  2.83  flow  over  a  20 
deg  compression  corner.  The  distribution  of  the 
peak  of  the  Reynolds  stress  profile  Is  exhibited  In 
Fig.  22.  It  shows  that  while  the  experimental  peak 
remains  approximately  constant  the  computed  peak 
steadily  diminishes.  In  Fig.  23,  the  corresponding 
distribution  of  the  location  of  the  peak  Ypea|<  of 
the  Reynolds  stress  profile  Is  displayed.  A 
pronounced  underprediction  of  the  distance  of  the 
peak  from  the  wall  Is  evident,  similar  to  that 
observed  In  Fig.  18  for  the  16  deg  corner. 

Discussion  of  Comparison  of  Reynolds  Stress 
Profiles -  - - 

In  summarizing  the  above  comparison  of 
Reynolds  shear  stress,  the  principal  discrepancy  Is 
the  underprediction  of  the  height  of  the  peak  of 
the  Reynolds  shear  stress.  The  Baldwln-Lomax  model 
is  modestly  successful  In  predicting  the  magnitude 
of  the  peak  of  the  Reynolds  shear  stress,  although 
the  success  Is  tempered  for  a  ■  20  deg  by  an 
apparent  Incorrect  trend  In  X. 

Recognizing  the  Inherent  simplicity  and 
limitations  of  the  mixing  length  concept.  It  is 
interesting,  nonetheless,  to  attempt  to  treat  the 
defects  of  the  8aldw1n-Lomax  model  "symptoma¬ 
tically".  It  is  noted  In  Fig.  23  that  the  height 
of  the  peak  Tpeak  of  the  computed  Reynolds  shear 
stress  correlates  with  the  magnitude  of  the 
computed  outer  length  scale  Ymax  of  the 
Baldwln-Lomax  model  for  a  *  20  deg;  a  similar 
observation  applies  for  a  ■  16  deg.  The  location 
of  the  experimental  peak  Reynolds  stress 
corresponds  to  the  outer  portion  of  the  boundary 
layer  ( 1 . e . ,  outside  the  point  where  the 
Baldwln-Lomax  model  switches  from  the  Inner  to  the 
outer  formulation).  This  suggests,  therefore,  that 
the  computed  Reynolds  shear  stress  may  be  improved 
by  Increasing  Ymax.  This  approach  was  attempted 
by  Vlsbal^S  for  a  *  16  deg.  Specifically,  was 
kept  constant  at  its  upstream  value;  this 
represents  an  Increase  In  Ymax  compared  to  the 
calculations  with  the  relaxation  eddy  viscosity 
model  (Fig.  24)  The  effect  of  Increasing  Ymaj(  is 
seen  In  Figs.  17  and  18.  The  magnitude  of  the  Deak 
Reynolds  shear  stress  Is  overpredlcted  for  x  <  Zim, 
while  a  slight  improvement  in  Yoeai<  is  noted.  The 
computed  Reynolds  shear  stress  profiles4®,  display 


a  double-peaked  oehavior  at  t/lm  *  0.49  and  0.98, 
in  disagreement  with  experiment.  It  may  be 
concluded,  therefore,  that  no  overal  1  improvement 
was  obtained  by  Increasing  Ymax. 

It  Is  evident  that  the  simple  mixing-length 
Baldwln-Lomax  model  is  incapable  of  accurately 
predicting  the  reattachment  and  downstream  recovery 
of  a  separated  2-D  compression  corner  flow.  The 
model  Is  based  upon  the  concept  of  an  equilibrium 
turbulent  boundary  layer  exhibiting  one 
characteristic  velocity  scale4^.  Qownstream  of 
reattachment,.  there  are  two  characteristic  velocity 
scales  of  the  turbulence,  namely  1)  an  outer 
velocity  scale  associated  with  the  turbulence 
fluctuations  In  the  outer  portion  of  the 
reattaching  free  shear  layer,  ana  2)  an  inner 
velocity  scale  u»  =  {  **(  x  ) /ow(  *  )  )  "w  associated 
with  the  Imposition  of  the  no  slip  bounaary 
condition  downstream  of  reattachment,  which  creates 
an  "inside  layer"  within  the  boundary  layer.  The 
failure  of  the  "simple"  extension  to  the 
Baldwln-Lomax  model  described  above  is  therefore 
not  surprising  within  this  framework.  A  more 
physical  realistic  turbulence  model  Is  required 
for  2-0  separated  compression  corner  flows  wnich 
Incorporates  the  effect  of  the  upstream  history  on 
the  turbulent  flow  and  the  oscillatory  motion  of 
the  shock  wave  structure43,44  The  assumption  of 
an  algebraic  eddy  viscosity  model  precludes  the 
Incorporation  of  the  turbulence  history,  except 
through  the  crude  technique  of  the  relaxation 
model.  With  regard  to  the  unsteady  large  amplitude 
Shock  wave  motion,  additional  research  is  needed  to 
elucidate  Its  effect  on  the  turbulence  structure 
and  on  the  recovery  of  the  boundary  layer 
downstream  of  reattachment. 


IV.  Conclusions 

A  comparative  study  has  been  performed  for  the 
MacCormack  hybrid  and  the  Beam-Warming  fully 
Implicit  algorithms  for  a  shock  wave- 1 urbu  1  en t 
boundary  layer  interaction  over  a  two-dimensional 
corner.  The  computations  employed  identical  grids 
and  the  Baldwln-Lomax  algebraic  turbulent  eddy 
viscosity  model.  It  is  observed  that  the  steady 
state  solution  of  the  MacCormack  hybrid  algorithm 
Is  remarkably  insensitive  to  Courant  number.  The 
accuracy  of  the  steady  state  solution  using 
MacCormack's  hybrid  algorithm  is  comparable  to  that 
of  the  Beam-Warming  method  for  all  cases.  Based 
on  experience  with  the  Mach  1.96  computations.  The 
MacCormack  hybrid  method  Is  observed  to  reduce  the 
computing  time  by  a  factor  of  up  to  3  relative  to 
the  Beam-Warming  method. 

The  computed  Reynolds  stress  profiles  are 
compared  with  the  experimental  data  of  Muck  et 
a]30-32.  it  Is  noted  that  the  magnitude  of  the 
peak  of  the  computed  Reynolds  shear  stress  is  in 
approximate  agreement  with  the  measurements, 
although  an  apparent  Incorrect  trend  is  evident  for 
a  =  20  deg.  The  major  discrepancy  is  the 
underprediction  of  the  location  of  the  peak  of  the 
computed  Reynolds  shear  stress.  It  is  noted  that  a 
simple  modification  of  the  Baldwln-Lomax  turbulence 
model  involving  an  increase  in  the  length  scale 
ymax  of  The  outer  eddy  viscosity  fails  to 
demonstrate  overall  improvement.  The  Baldwln-Lomax 
model  is  based  on  the  mi x 1 ng- length  concept,  and  is 
incapable  of  accurately  predicting  the  recovery  of 
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a  separated  2-0  compression  corner  flow.  It  Is 
noted  that  several  additional  physical  factors, 
omitted  from  the  theoretical  model,  also  affect  the 
recovery  of  the  boundary  layer  Including  the 
history  effect  of  the  turbulence  structure  and  the 
large  amplitude  oscillatory  motion  of  the  shock 
structure. 
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Fig.  2  -  Surface  pressure  for  Mach  1.96 
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Fig.  3  -  Skin  friction  for  Mach  1.96  flow 
over  16  deg  ramp 
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rig.  16  -  Reynolds  shear  stress  at  X/ >  =  5.4 
for  Mach  2.83  flow  over  16  deg  ramp. 


fig.  19  -  Reynolds  shear  stress  at  X/;»  =  1.0 
for  Mach  2.83  flow  over  20  deg 
ramp. 


- CFL  =  72. 

- CFL  =  30. 

---  BERM- WORMING 
□  EXPERIMENT 


~-4  .  / 

5 

4.0  "  Q 


I.  2.0 


j 

- CFL  s  85. 

-  -  CFL  s  30.  ^ 

- 9ERM-WORM  1RELOXOTI0N1 

- 8EOM-WORM  I CONST  TMOX1 

a  EXPERIMENT 


0.0  *— 
0.0 


0.0  <-* 

-3.0 


x  /  S, 


Fig.  17  -  Magnitude  of  peak  of  the  Reynolds 

stress  for  Mach  2.83  flow  over  16  deg 
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Fig.  20  -  Reynolds  shear  stress  at  X/ '>  =  2.1 
for  Mach  2.83  flow  over  20  dea  ramo. 
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Fig.  18  -  Location  of  peak  in  Reynolds  shear 
stress  profile  for  Mach  2.83  flow 
over  16  deg  ramp. 
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Fig.  21  -  Reynolds  shear  stress  at  2  4.6 

for  Mach  2.83  flow  over  20  -jeo 
ramp . 
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Fig.  22  -  Magnitude  of  peak  of  the  Reynolds 
stress  for  Mach  2.83  flow  over 
20  deg  ramp. 
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Fig.  24  -  Baldwin-Lomax  length  scale  Ymax  for 
Mach  2.83  flow  over  16  deg  ramp. 
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Fig.  23  -  Baldwin-Lomax  length  scale  Ymax, 
and  location  of  the  peak  of  the 
Reynolds  shear  stress  for  Mach  2.83 
flow  over  20  deg  ramp. 
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Abstract 

The  3-0  shock  wave-turbulent  boundary  layer 
interaction  generated  by  a  sharp  fin  Is  examined 
both  experimentally  and  theoretically  at  Mach  3  for 
a  fin  angle  ag  *  20  deg  and  Reynolds  number  Reg  * 
9  x  105.  This  study  represents  an  extens1on"of 
previous  research  for  the  sharp  fin  configuration 
to  stronger  Interactions.  The  experimental  data 
Include  surface  pressure  profiles,  surface 
streamline  patterns,  and  boundary  layer  profiles  of 
pitot  pressure  and  yaw  angle.  Two  separate 
theoretical  approaches  or  "models"  were  employed. 
Both  models  employ  the  3-0  compressible  Navler- 
Stokes  equations  In  mass-averaged  variables.  The 
theoretical  approach  of  Knight  utilizes  the 
algebraic  turbulent  eddy  viscosity  model  of  Baldwin 
and  Lomax,  and  the  theoretical  model  of  Horstman 
employs  the  two-equation  turbulence  model  of  Jones 
and  Launder  coupled  with  the  wall  function  model  of 
Vlegas  and  Rubesln.  The  coeluted  surface 
pressure,  surface  streamlines,  pitot  pressure  and 
yaw  angle  profiles  are  in  good  agreement  with  the 
experimental  data,  thereby  confirming  the  efficacy 
of  the  theoretical  approaches  which  were  previously 
validated  for  the  3-0  sharp  fin  conf Iguratlon  at 
Mach  3  for  smaller  ag  ( 1 . e. ,  weaker  Interactions). 
The  three  dimensional  velocity  fields  computed  by 
both  models  are  In  close  agreement,  although  the 
eddy  viscosity  profiles  differ  significantly  within 
the  3-0  Interaction.  This  result  Indicates  that 
the  overall  structure  of  this  3-0  sharp  fin 
Interaction  Is  Insensitive  to  the  turbulence  model. 
A  series  of  particle  pathllne  traces  were  examined 
for  each  model,  and  found  to  be  In  close  agreement. 
The  calculated  flowflelds  display  a  prominent 
vortical  structure  associated  with  the  shock¬ 
boundary  layer  Interaction  In  agreement  with  the 
flowfleld  models  of  Token,  and  Kubota  and  Stollery. 
The  structure  Is  characterized  by  two  significant 
surfaces,  namely,  1)  the  surface  of  separation 
which  emanates  from  the  line  of  coalescence  (line 
of  separation)  and  spirals  Into  the  vortical  core, 
and  2)  the  surface  of  attachment  which  intersects 
the  wall  at  the  line  of  divergence  (line  of 
attachment)  and  demarcates  the  fluid  which  Is 
entrained  Into  the  vortical  structure. 
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I .  Introduction 

For  nearly  forty  years,  significant  Interest 
has  been  focused  on  the  subject  of  two-  and 
three-dimensional  shock  wave-turbulent  boundary 
layer  Interactions  (denoted  “2-0"  or  "3-0  turbulent 
Interactions"  for  brevity).  The  phenomenon  Is 
widespread  In  aerodynamics,  turbomachinery  and 
other  areas  of  fluid  mechanics,  and  therefore  a 
clear  motivation  exists  for  developing  a  deeper 
understanding.  A  variety  of  simplified  geometrical 
configurations  have  been  employed  to  Investigate 
3-0  turbulent  Interactions*.  A  sample  of  these 
Include,  1)  sharp  fin  mounted  perpendicular  to  a 
flat  plate,  2)  blunt  fin  mounted  perpendicular  to  a 
flat  plate,  3)  swept  compression  corner,  4)sem1cone 
affixed  to  a  flat  plate,  5)  cone  within  a  circular 
wind  tunnel,  and  6)  normal  wall  Jet. 

The  focus  of  the  present  paper  Is  the  3-0 
oblique  shock  wave- turbulent  boundary  layer  Inter¬ 
action  generated  by  a  sharp  fin  attached  to  a  flat 
plate  (Fig.  1).  The  overall  flowfleld  Is  determined 
by  a  small  set  of  parameters,  namely, 
the  upstream  Mach  number  M»,  the  Reynolds  number 
Re(«  based  upon  the  boundary  layer  thickness  at 
the  strearawlse  station  corresponding  to  the  leading 
edge  of  the  fin  (where  im  Is  measured  in  the 
absence  of  the  fin),  the  nature  of  the  thermal 
boundary  condition  on  the  flat  plate  and  fin  (e.g., 
adiabatic  or  fixed  temperature),  and  the  fin  angle 
ag.  This  configuration  has  been  the  subject  of 
several  experimental  and  theoretical  Investi¬ 
gations.  Experiments  have  focused  principally  on 
surface  measurements,  and  Include  the  studies  of 
Stanbrook2,  McCabe3,  Law4,  Kubota  and  Stollery5, 
Zheltovodov3,  Oolllng?  and  Goodwin3.  In  recent  years, 
detailed  boundary  layer  measurements  have  been 
obtained  by  Peake",  Oskam,  Vas  and  Bogdonoff *0_*2 , 
and  McClure  and  Dolling*3.  Numerical  simulations 

using  the  3-0  Reynolds-averaged  Navler-Stokes 

equations  have  been  performed  principally  at  Mach  3, 
and  Include  the  studies  of  Horstman  and  Hung*4  and 
Knight*5"*'.  These  computations,  detailed  In  Table 
1,  have  previously  considered  fin  angles  an  up 
to  10  deg.  The  Investigations  of  Horstman 
and  Hung  utilized  the  Escudler*3  turbulence  model, 
while  the  later  work  of  Horstman*3  employed  the 
Jones-Launder^O  model.  The  calculations  of  Knight 
utilize  the  Baldwln-Lomax^l  turbulence  model. 

These  prior  calculations  have  been  examined  in 
comparison  with  experimental  data  for  a  wide 
variety  of  flow  quantities  Including  surface 


pressure,  skin  friction,  heat  transfer,  pitot 
pressure,  yaw  angle,  pitch  angle  and  static 
pressure.  In  general,  good  agreement  was  obtained 
with  the  experimental  data. 


II.  Description  of  Experiment 

The  experiments  were  performed  In  the 
supersonic  high  Reynolds  number  wind  tunnel  at  the 


TABLE  1. 


Computations  of  3-0  Sharp  Fin  Configuration  at  Mach  3 


4 


Investigator 

ftOi  P2/P«  Theoretical  Experimental 

(d«J) 


3.75 

8.7 

X 

10f 

1.32 

Horstman  and  Hung14 

Oskam  et  a  1. 10-12 

3.75 

9.3 

X 

105 

1.32 

Knight15 

Oskam  et  a l.10-12 

10 

2.8 

X 

105 

2.01 

Knight16 

McClure  and  Oolling^ 

10 

3.4 

X 

lof 

2.01 

Horstman1  ^ 

McClure  and  Dolllnq12 

10 

8.7 

X 

105 

2.01 

Horstman  and  Hung14 

Oskam  et  al . 10-12 

10 

9.3 

X 

105 

2.01 

Knight15 

Oskam  et  al. 10-12 

20 

8.6 

X 

105 

3.65 

Knight* 

Shapey  and  Bogdonoff*. 
Goodwl n® 

20 

3.6 

X 

105 

3.65 

Horstman* 

Shapey  and  Bogdonoff*. 
Goodwin® 

'Present  paper 

Note:  Actual  Mach  number  Is  2.94  t  .01. 


There  are  three  objectives  for  the  present 
paper: 

1.  Examine  the  Accuracy  of  the  Theoretical 
Models  for  Stronger  Interactions. 

As  indicated  in  Table  1,  the  strongest  3-0 
sharp  fin  Interaction  computed  previously  at  Mach  3 
corresponded  to  a  pressure  ratio  P2/P»  *  2.01  (o^  * 
10  deg),  where  p.  Is  the  upstream  static  pressure 
•  and  P2  Is  the  theoretical  downstream  Invlscld 
pressure.  A  critical  Issue  Is  the  examination  of 
the  accuracy  of  the  theoretical  models  for  stronger 
Interactions.  In  the  present  study,  the 
theoretical  models  are  examined  for  the  3-0  sharp 
fin  Interaction  at  Mach  3  and  og  *  20  deg,  which 
exhibits  a  pressure  rise  of  3./.  This  is  the 
strongest  Interaction  considered  for  these  models 
at  Mach  3. 


Princeton  University  Gas  Dynamics  Laboratory.  The 
facility  has  a  20  cm  x  20  cm  test  section,  with  a 
nominal  freestream  Mach  number  of  2'.  93.  The 
settling  chamber  pressure  and  temperature  were  6.8 
x  IQ5  Pa  tit  and  251  *K  t5t,  respectively,  yielding 
a  nominal  Reynolds  number  of  7.0  x  10®  m-1.  The 
experiments  were  performed  under  near  adiabatic 
wall  conditions. 

The  sharp  fin  Is  14.21  cm  long  and  12.7  cm 
high.  The  fin  was  fabricated  from  alumlmxn  with  a 
Sharp  unswept  leading  edge,  and  oriented  at  a  right 
angle  to  the  tunnel  wall  ("flat  plate”).  The  fin 
was  mounted  In  a  unique  variable-geometry  apparatus 
which  permitted  the  achievement  of  fin  angles 
exceeding  20  degrees,  thereby  extending  the  range 
of  the  experiments  beyond  the  earlier  fixed- 
geometry  configuration10-12. 


2.  Comparison  of  two  different  Theoretical 
Models. 

An  important  objective  of  the  present  research 
Is  to  examine  the  computed  flowflelds  for  the  3-0 
sharp  fin  Interaction  obtained  using  two  different 
turbulence  models.  In  this  effort,  Navler-Stokas 
calculations  have  been  performed  by  Knight  and 
Horstman  using  the  Baldwln-lomax  and  Jones-launder 
turbulence  models,  respectively,  for  the  3-0  sharp 
fin  configuration.  A  principal  Issue  Is  the 
determination  of  the  sensitivity  of  the  computed 
flowfleld  to  the  turbulence  model  employed. 

3.  Examine  the  Flowfleld  Structure  of  the  3-0 
Sharp  Fin  Interaction. 

Provided  the  theoretical  models  yield  good 
agreement  with  the  experimental  data  for  the  Mach 
3,  a-  *  20  deg  configuration,  the  computed 
flowflelds  can  be  utilized  to  examine  and 
understand  the  flowfleld  structure  of  this  3-0 
sharp  fin  Interaction.  Flowfleld  models  have  been 
developed,  for  example,  by  Token22,  and  Kubota  and 
Stollery5.  These  models  may  be  examined  using  the 
computed  flowflelds. 


Surface  pressure  distributions  were  obtained 
along  rows  of  orifices  aligned  with,  the 
x-dlrectlon.  A  kerosene- lampblack  technique21  was 
employed  to  obtain  surface  flow  angularity.  The 
boundary  layer  on  the  tunnel  wall  ("flat  plate”) 
was  surveyed  using  a  computer-controlled  nulling 
cobra  probe12  which  measured  pitot  pressure  p0  and  * 
yaw  angle  s,  where  a  *  tan-1(w/u)  with  (u,v,w) 
Indicating  the  cartesian  velocity  components  In  the 
( x ,y ,z )  coordinate  system  (Fig.  1).  The  survey 
locations,  shown  In  Fig.  2,  were  selected  to 
provide  detailed  Information  within  the  region 
between  the  line  of  coalescence  ("3-0  separation 
line"’)  and  shock  wave,  and  the  region  downstream 
of  the  shock.  The  specific  survey  stations  are 
Indicated  In  Table  2,  where  xs  »  x  -xsn0ck  (*shock 
Is  the  streamwlse  location  of  the  theoretical 
Invlscld  shock  wave  at  the  specified  soanwlse 
location  z),  z  Is  the  spanwlse  location,  and  4-  * 
1.27  cm. 

The  Incoming  flow  on  the  tunnel  wall  is  an 
equilibrium,  two-dimensional  boundary  layer  which 
has  been  extensively  surveyed24*25  amj  observed  to 
closely  fit  the  law  of  the  Wall  and  Wake26.  The 
thickness  of  the  Incoming  boundary  layer  4.  Is 
approximately  1.4  cm. 
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TABLE  2.  Boundary  Layer  Survey  Locations 


Survey  Location 

xs/«« 

Z/4. 

1 

-5.40 

5.81 

2 

-4.40 

5.81 

3 

-3.40 

5.81 

4 

-2.40 

5.81 

5 

-1.40 

5.81 

6 

-0.40 

5.81 

7 

0.60 

5.81 

a 

2.60 

5.81 

9 

-3.94 

7.81 

10 

-0.14 

4.81 

li 

1.13 

3.81 

III.  Description  of  Computations 
A.  Theoretical  Approach  of  Knight 

1.  Governing  Equations  and  Numerical  Algorithm 

The  governing  equations  are  the  full  mean 
compressible  3-0  Navler-Stokes  equations  using 
mass -averaged  variables27  and  strong  conservation 
form28.  The  molecular  dynamic  viscosity  Is  given 
by  Sutherland's  law.  The  molecular  and  turbulent 
Prandtl  numbers  are  0.73  and  0.9,  respectively. 
Turbulence  Is  modelled  through  the  inclusion  of  the 
two-layer  algebraic  turbulent  eddy  viscosity  model 
of  Baldwin  and  Lomax23,  with  the  mixing  length 
specified  by  the  formula  of  Bulaev29  as  discussed 
In  Gessner  and  Po30.  The  Baldwln-Lamax  (BL)  model 
Is  Implemented  as  discussed  In  Ref.  15-17. 

A  3-0  coordinate  transformation  (C(x,y,z), 
n(x,y,z),  c(x,y,z))  Is  used  to  map  the  physical 
domain  (shown  as  the  dotted  lines  In  Pig.  1)  Into  a 
cube  In  the  transformed  domain  whose  simple  shape 
facilitates  the  application  of  the  numerical  algo¬ 
rithm.  Various  methods  have  been  developed  for 
numerical  generation  of  curvilinear  coordinates33. 
For  the  3-0  sharp  fin,  however,  a  simple  analytic 
transformation  was  employed  (see  below). 

At  the  upstream  boundary  (A8H6  In  Fig.  1),  the 
flow  variables  are  held  fixed  at  the  values  corres¬ 
ponding  to  a  developed  boundary  layer  whose  proper¬ 
ties  are  In  close  agreement  with  the  experiment.  On 
the  solid  surfaces  (A8C0EF  and  FEKL)  the  velocity  Is 
zero,  a  fixed  surface  temperature  (near  adiabatic 
conditions)  is  applied,  and  the  normal  derivative  of 
the  static  pressure  Is  set  to  zero.  On  the  plane  of 
symmetry  (AFLS) ,  the  normal  component  of  the  veloc¬ 
ity  Is  set  to  zero,  and  the  normal  derivatives  of 
the  remaining  flow  quantities  are  zero.  The  right 
boundary  BCDJIH  is  sufficiently  far  from  the  fin  to 
Insure  that  the  boundary  layer  Is  locally  2-0,  and 
therefore  a  simple  gradient  boundary  condition  »/»z 
■  0  Is  employed.  At  the  downstream  boundary,  the 
conventional  i/H  •  0  condition  Is  specified. 

The  governing  equations  are  solved  by  an 
efficient  hybrid  explicit-implicit  numerical 
algorithm3®.  The  technique  utilizes  the  second- 
order  accurate  explicit  method  of  MacCormack32*33, 
and  the  second-order  accurate  Implicit  method  ("Box 
Scheme")  of  Keller34.  The  Implicit  algorithm  of 
Keller  Is  employed  In  a  thin  layer  (denoted  the 
"computational  sublayer")  adjacent  to  the  solid 


boundaries  where  the  large  flow  gradients  require 
exceptionally  fine  grid  spacing  for  accurate 
resolution.  The  Box  Scheme  Is  applied  to  the 
asymptotic  form  of  the  Navler-Stokes  equations  in 
this  region,  whose  height  is  restricted3®-37*3®-37 
according  to  the  expression  s  60,  where  z^+ 
*  ZinU»/vw,  Is  the  local  height  of  the 
computational  sublayer.  u«  Is  the  local  friction 
velocity  (u*  »  (s^/ow)^)*  and  vw  1s  the  kinematic 
molecular  viscosity  evaluated  at  the  surface.  This 
layer  Is  typically  a  few  percent  of  the  local 
boundary  layer  thickness.  The  explicit  algorithm 
of  NacCormack  Is  applied  to  the  full  Navler-Stokes 
equations  In  the  remainder  of  the  flowfleld 
(denoted  the  "ordinary  region"). 

The  hybrid  algorithm  has  been  successfully 
applied  to  a  wide  range  of  two-  and  three- 
dimensional  flows  exhibiting  shock-boundary  layer 
interaction  and  flow  separat1onl5-  17,35-37.  rhe 
code  Is  written  In  CYBER  200  FORTRAN,  and  executes 
on  the  VPS  32  at  NASA  Langley  Research  Center.  The 
VPS  32  Is  a  vector-processing  supercomputer  which 
Is  archl  tectural  ly  similar  to  the  CYBER  205.  The 
explicit  portion  of  the  algorithm  is  highly 
vectorized,  with  typical  vector  lengths  of  1500, 
and  has  achieved  an  execution  rate  of  approximately 
100  MFlops  (million  floating  point  operations  per 
second)  on  the  VPS  32  using  a  32-bit  word  length. 

2.  Details  of  Computation 

The  upstream  boundary  layer  profile  was 
computed  using  a  separate  boundary  layer  code38. 
The  flow  conditions  are  Indicated  in  Table  3,  and 
are  In  close  agreement  with  the  experimental 
conditions  and  upstream  profile  of  Horstman. 

The  finite-difference  mesh  was  generated 
according  to  the  method  described  in  Ref.  15.  A 
total  of  32  streamwlse  grid  planes  were  utilized, 
uniformly  spaced  In  the  x-dlrectlon  with  ax  3  4.. 
The  upstream  boundary  was  located  at  56*  upstream 
of  the  fin  leading  edge,  and  the  downstream 
boundary  at  x  *  26a*.  The  grid  spacing  within  each 
plane  was  a  combination  of  geometrically-stretched 
and  uniformly  spaced  points.  The  number  of 
ordinary  points  In  the  y-and  z-dlrectlons  are  32 
and  48,  respectively.  The  computational  sublayer 
was  resolved  using  8  points  in  the  direction  normal 
to  the  surface.  A  separate  refined  grid  was 
utilized  In  the  sublayer  region  in  the  inuiedlate 
neighborhood  of  the  corner  formed  by  the  flat  plate 
and  the  fin.  The  total  number  of  grid  points  was 
64,956.  The  height  of  the  first  grid  point 
adjacent  to  the  fin  or  flat  plate  was  less  than  3.0 
wall  units  at  all  location  (l.e.,  az£*  s  3.0, 
where  az£+  »  az£u*/vw'  and  az£  is  the 
distance  of  the  first  row  of  grid  points  adjacent 
to  the  surface).  Two  separate  computations  were 
performed  to  examine  the  sensitivity  of  the 
solution  to  the  height  z^  of  the  computational 
sublayer  adjacent  to  the  flat  plate.  These 
computations  employed  z,),  »  5.33  x  10*3  cm  and 
7.62  x  10"3  cm.  The  maximum  grid  spacing  In  the 
y-dlrectlon  for  these  cases  was  Ay.  3  0.58*.  and 
0.59a«.  The  height  of  the  computational  domain  was 
8a..  The  width  of  the  domain  Increased  linearly 
from  13.04.  at  x  3  0  to  32.64.  at  x  3 26a..  The 
maximum  grid  spacing  In  the  z-dlrectlon  varied 
between  0.424.  and  1.07a..  The  comouted  results 
using  the  two  separate  grids  were  found  to  be 
essentially  Identical. 
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TABLE  3.  Flow  Conditions  for  3-0  Sharp  Fin  at  ag  •  20  deg 


Case 

<• 

(cm) 

pt. 

(kP») 

(deg*K) 

Experiment  (Shapey) 

1.4 

2.93 

9.8  x  105 

690 

251 

Theory  (Knight) 

1.3 

2.93 

8.8  x  IQS 

690 

256 

Theory  (Horstnan) 

1.4 

2.94 

8.8  x  IQ* 

690 

267 

B.  Theoretical  Modal  of  Horstnan 

1.  Governing  Equations  and  Numerical  Algorithm 

The  governing  equations  are  the  full  mean 
compressible  3-0  Navler-Stokes  equations  using 
mass-averaged  variables.  The  molecular  dynamic 
viscosity  Is  specified  using  Sutherland's  law.  The 
molecuiar  and  turbulent  Prandtl  nixnbers  are  0.72 
and  0.90,  respectively.  The  effects  of  turbulence 
are  modelled  using  the  two-equation  eddy  viscosity 
model  of  Jones  and  launder^0  (JL).  An  additional 
element  of  the  turbulence  modelling  Is  the 
Incorporation  of  the  compressible  wall  functions  of 
Vlegas  and  Rubes1n39,  previously  utilized  In  the 
study  of  the  three  dimensional  shock  boundary  layer 
Interactions  for  a  swept  compression  corner 
configuration40.  The  other  boundary  conditions 
employed  are  similar  to  those  described  previously. 

The  governing  equations  are  solved  by  the 
explicit  numerical  algorithm  of  NacCormack32.  The 
algorithm  has  been  widely  employed  for  the 
computation  of  2-0  and  3-0  turbulent  Interactions. 
The  nunerlcal  code  was  executed  on  the  CRAY  X-MP/22 
at  NASA  AMES.  The  code  Is  fully  vectorized  using 
CRAY  FORTRAN  and  utilizes  the  CRAY  Solid  State 
Disk. 

2.  Petal  Is  of  Computations 

The  upstream  boundary  layer  profile  was 
computed  using  a  separate  boundary  layer  code 
employing  the  two-equation  W1 1 cox-Rubes In 
turbulence  model4!.  The  upstream  flow  conditions 
are  indicated  In  Table  3,  and  are  In  close 
agreement  with  the  experimental  conditions  and  the 
upstream  profile  of  Knight. 

The  numerical  grid  was  generated  algebraically 
using  a  combination  of  geometric  stretching  and 
uniform  spacing.  A  total  of  64  streamwlse  grid 
planes  were  employed,  spaced  uniformly  with  ax  « 
0.39a..  The  upstream  boundary  was  located  1.8(. 
upstream  of  the  fin  leading  edge,  and  the 
downstream  boundary  was  positioned  it  *  *  21.94.. 
Within  each  streamwlse  grid  plane,  the  total  masher 
of  grid  points  In  the  y-  and  z-dlrectlons  was  32 
and  44,  respectively.  The  maximum  grid  spacing  In 
the  y-  and  z-dlrectlons  was  0.39l«.  The  height  of 
the  computational  domain  In  the  y-dlrectlon  was 
7.2i«  and  the  width  (measured  from  the  plane  of 
syemetry  or  the  fin  surface)  was  13.04..  The  total 
number  of  grid  points  was  90,112. 

IV.  Results 

A.  Comparison  with  Experiment 

The  computed  and  measured  pitot  pressure 
profiles  at  Station  1  (not  shown),  located  upstream 


of  the  Interaction,  are  In  close  agreement. 
Similarly,  the  computed  and  experimental  pitot 
profiles  at  Station  2  (not  shown),  located  at  the 
line  of  upstream  Influence  (as  defined  by  the 
experimental  surface  pressure),  are  in  close 
agreement,  and  exhibit  negligible  deviation  from  an 
equilibrium  2-0  profile.  The  comouted  and 
experimental  pitot  pressure  profiles  at  Stations  3 
through  8,  10  and  11  are  displayed  In  Figs.  3  to 
10.  The  horizontal  axis  Is  the  pitot  pressure  p0, 
normalized  by  the  upstream  freestream  pitot 
pressure  Pp^.  The  vertical  axis  is  the  distance 
measured  from  the  flat  plate,  normalized  by  the 
upstream  boundary  layer  thickness  4..  It  is  noted 
that  the  upstream  boundary  layer  thickness  is  not. 
In  general,  the  appropriate  vertical  scaling 
parameter  for  this  interaction.  The  experimental 
data  of  McClure  and  0ollingl3  suggest  that  the 
appropriate  vertical  scaling  is  given  by 

yRe^/’/40,  where  40  Is  the  experimental 
boundary  layer  thickness  measured  immediately 
upstream  of  the  shock  (with  the  fin  removed)and  at 
the  specified  spanwlse  location.  The  choice  of  4. 
as  the  vertical  scaling- was  motivated  by  the  desire 
to  clearly  portray  the-  vertical  extent  of  the 
Interaction  relative  to  the  height  of  the  upstream 
boundary  layer.  Similarly,  the  profiles  are  shown 
at  selected  values  of  xs/4«,  where  xs  »  x  -  xshock 
and  x$hock  *s  the  location  of  the  theoretical 
Invlscld  shock  at  the  specified  spanwise  location. 

It  Is  noted  that  the  observations  of  Settles  and 
Bogdonoff4^,  Dolling  and  Bogdonoff4!,  and  Lu  and 
Settles44  Indicate  that  the  appropriate  scaling  is 
given  by  xj  *  *sftelo3/<o- 

In  Fig.  3,  profiles  of  pitot  pressure  are 
shown  at  Station  3,  which  is  coincident  with  the 
line  of  coalescence  as  defined  by  the  kerosene 
lampblack  visualization.  Although  displaying  a 
slightly  greater  pitot  pressure  In  the  outer 
portion  of  the  boundary  layer,  the  computed  results  * 
are  In  reasonable  agreement  with  experiment. 
Similar  agreement  Is  obtained  at  Station  9  (not 
shown).  In  Fig.  4,  the  results  are  displayed  at 
Station  4,  located  approximately  one-third  of  the 
distance  between  the  line  of  coalescence  and  the 
shock  wave  at  this  spanwlse  position.  The  computed 
and  experimental  profiles  display  a  modest 
"overshoot"  outside  the  boundary  layer,  associated 
with  the  compression  system  ahead  of  the  shock 
wave!3,I6,17.  jn  Fig.  5,  the  results  are  displayed 
at  Station  S,  located  approximately  two-thirds  of 
the  distance  between  the  line  of  coalescence  and 
the  theoretical  Invlscld  shock  at  this  spanwise 
position.  The  experimental  profile  displays  a 
slight  S-shaped  behavior  near  the  wall,  which  is 
less  apparent  In  the  computed  profiles.  The 
overshoot  In  pitot  pressure  Is  more  pronounced  at 
this  location.  In  both  figures,  the  computed 
profiles  are  In  reasonable  agreement  with  the 
experimental  data,  and  accurately  predict  the 


observed  overshoot  In  pp.  The  shock  capturing 
nature  of  both  menerlcal  algorithms  can  be  seen  In 
the  smearing  of  the  pitot  pressure  profile  near 
y  *  1-5*.. 

The  calculated  and  experimental  pitot  pressure 
profiles  at  Station  6  are  shown  In  Fig.  6.  Que  to 
the  close  proximity  of  this  station  to  the  shock 
wave,  uncertainty  exists  In  the  measurement  of 
pitot  pressure  and  yaw  angle  outside  the  boundary 
layer,  and  the  experimental  data  has  therefore  been 
denoted  by  a  dotted  line  for  y  >  1.54®.  Further 
experimental  Investigation  Is  required  to  resolve 
this  Issue.  Within  the  boundary  layer,  reasonable 
agreement  Is  obtained  between  the  computation  and 
measurement.  The  S-shape  character  of  the  profile 
Is  again  apparent,  with  reasonable  agreement 
between  theory  and  experiment.  The  region  of  high 
pitot  pressure  near  the  surface  Is  associated  with 
a  local  maximum  In  the  Mach  number.  In  Fig.  7, 
pitot  profiles  are  displayed  at  Station  10,  located 
inmedlately  upstream  of  the  shock  and  closer  to  the 
fin.  The  experimental  data  outside  the  boundary 
layer  Is  again  subject  to  uncertainty  due  to  the 
proximity  of  the  shock  wave.  Within  the  boundary 
layer,  the  agreement  between  the  theory  and 
experiment  Is  good. 

Pitot  pressure  results  at  Stations  7  and  11, 
located  Immediately  downstream  of  the  shock  wave, 
are  displayed  In  Figs.  8  and  9.  Good  agreement  Is 
again  observed.  The  discrepancy  In  the  computed 
pitot  pressure  outside  the  boundary  layer  In  Fig.  3 
Is  associated  with  the  shock-capturing  nature  of 
the  numerical  algorithms,  and  the  difference  In 
streamwlse  grid  spacing  for  the  two  computations, 
and  the  proximity  of  Station  7  to  the  shock  (xs  * 
0.60a,).  In  Fig.  10,  results  are  shown  at  Station 
8,  located  furthest  downstream  of  the  shock.  The 
comparison  between  computed  and  experimental 
results  Is  good. 

The  computed  and  experimental  yaw  angle 
profiles  at  Stations  1  and  2  (not  shown)  display 
negligible  values  (<  3  deg).  The  calculated  and 
measured  yaw  angle  profiles  at  Stations  3  through 
8,  10  and  11  are  displayed  In  Figs.  11  to  18.  In 
Fig.  3,  the  yaw  profiles  are  shown  at  Station  3  (xs 
*  -3.4a,,  z  *  5 . 31  a,) ,  located  at  the  experimental 
line  of  coalescence.  It  Is  observed  that  the 
computed  profiles  near  the  surface  underpredict  the 
observed  yaw  angle.  This  Is  attributable  to  the 
differences  between  computed  and  experimental  lines 
of  coalescence.  In  particular,  the  computed  lines 
of  coalescence  for  the  8a ldw In-Lomax  (BL)  model  and 
Jones-Launder  (JL)  modal  are  located  at  xs  ■  -3.0a, 
and  -2.4s„  respectively,  at  z  *5.8la,.  At  Station 
3,  the  calculated  values  of  the  surface  yaw  angle 
are  37  deg  (BL)  and  12  deg  (JL),  while  the 
experimental  surface  yaw  angle  (based  on  kerosene 
lampblack  visualization)  Is  approximately  54  dog. 
A  similar  observation  applies  to  Station  9  (not 
shown). 

In  Fig.  12,  yaw  angle  profiles  are  shown- at 
Station  4,  located  approximately  one-third  of  the 
distance  between  the  line  of  coalescence  and  shock 
at  this  spanwlse  position.  The  computed  yaw  angles 
are  In  reasonable  agreement  with  experiment,  except 
In  the  Immediate  vicinity  of  the  surface,  where  the 
computed  profiles  disagree  by  15-20%  from  the 
experiment.  In  particular,  the  computed  surface 
yaw  angles  are  63  deg  (BL)  and  44  deg  (JL),  while 
the  experimental  value  Is  54  deg  based  upon  cobra 


probe  measurements  and  kerosene  lampblack 
visualization. 

Taw  angle  results  at  Station  5,  located 
approximately  two-thirds  of  the  distance  between 
the  line  of  coalescence  and  the  shock  wave  at  this 
spanwlse  position,  are  displayed  In  Fig.  13.  The 
agreement  between  the  theory  and  experiment  is 
good,  although  the  Jones-Launder  model  overpredicts 
the  yaw  angle  In  the  outer  portion  of  the  boundary 
layer.  The  computed  surface  yaw  angles  are  61  deg 
(BL)  and  53  deg  (JL),  In  reasonable  agreement  with 
the  experimental  value  of  60  deg. 

In  Figs.  14  and  15,  yaw  angle  profiles  are 
shown  at  Stations  6  and  10,  immediately  upstream  of 
the  shock.  As  discussed  previously,  uncertainties 
exist  in  the  experimental  data  outside  the  boundary 
layer  due  to  the  proximity  of  the  shock  wave,  ana 
the  data  are  consequently  Identified  by  a  dotted 
line  In  that  region.  Overall  good  agreement  is 
observed  between  the  calculated  and  experimental 
results  within  the  boundary  layer.  The  calculated 
surface  yaw  angles  at  Station  6  are  60  deg  (8L)  and 

57  deg  (JL),  In  reasonable  agreement  with  the 
experimental  value  of  64  deg.  At  station  10,  the 
computed  surface  values  are  61  deg  (8L)  and  58  deg 
(JL),  In  close  agreement  with  the  measured  value  of 

58  deg. 


The  calculated  and  experimental  yaw  angle 
profiles  at  Stations  7  and  11,  located  immediately 
downstream  of  the  shock,  are  shown  in  Figs.  16  and 
17.  The  calculated  results  are  again  observed  to 
be  In  close  agreement  with  experiment,  although 
displaying  a  somewhat  less  full  profile  near  the 
surface  at  Station  7.  The  predicted  surface  yaw 
angles  are  61  deg  (BL)  and  59  deg  (JL)  at  Station 
7,  In  close  agreement  with  the  experimental  value 
of  64  deg.  At  Station  11,  the  calculated  surface 
values  are  59  deg  (BL)  and  51  deg  (JL),  and  the 
measured  value  Is  55  deg.  In  Fig.  18,  yaw  angle 
profiles  are  shown  at  Station  8,  downstream  of  the 
shock.  The  calculated  and  experimental  profiles 
are  observed  to  be  In  excellent  agreement.  The 
predicted  surface  values  of  52  deg  (BL)  and  50  deg 
(JL)  are  In  close  agreement  with  the  experimental 
value  of  48  deg. 

The  calculated  surface  pressure  for  both 
models  has  been  compared  with  the  data  of  Goodwin8 
for  the  same  configuration.  The  models  accurately 
predict  the  extent  of  the  upstream  Influence,  and 
the  pressure  distribution  from  the  upstream* 
Influence  location  to  the  plateau  region. 
Downstream  of  the  plateau  region,  the  computed 
pressures  moderately  underpredict  the  data. 


Further  Comparison  of  Computed  Flowfleld 


A  detailed  comparison  of  the  computed 
flowflelds  of  knight  and  Horstman  was  performed  to 
determine  the  extent  of  similarity  of  the  two 
theoretical  approaches.  Profiles  of  the  computed 
x-component  velocity,  yaw  angle,  pitch  angle  r  = 
tan'I  ( v/Vu^-t-w? ) ,  pitot  pressure,  Mach  number  and 
turbulent  eddy  viscosity  were  examined  at  a 
selected  streamwlse  station  x  *  11a®  for  i  -if\n  = 
a®  to  10a®  In  Increments  of  a®,  where  Zf1n  is  the 
width  of  the  fin  at  a  given  x.  A  representatl ve 
sample  Is  displayed  In  Figs.  19  to  22, 
corresponding  to  x  *  11a®  and  z  -  Zfin  *  6a®.  The 
position  is  roughly  halfway  between  the  line  of 
coalescence  and  the  shock  wave  (as  measured  in  the 


streamwlse  direction)  at  this  spanwlse  location. 
The  x-ccmponent  velocity,  yaw  angle  and  pitch  angle 
profiles,  shown  In  Figs.  19,  20  and  21, 
respectively,  are  observed  to  be  In  very  close 
agreement,  differences  In  the  computed  yaw  angle 
occur  only  within  the  region  y  <  0.3<*.  Note  that 
It  Is  difficult  to  define  a  "local"  boundary  layer 
thickness  within  the  3-0  Interaction  region  due  to 
the  non-uniformity  of  the  Invlscld  flowi0»^«12. 
In  Fig.  20,  however,  the  eddy  viscosity  profiles 
Indicate  a  significant  difference;  In  particular, 
the  peak  values  of  the  turbulent  eddy  viscosity  c 
differ  by  a  factor  of  fourteen.  It  Is  emphasized 
that  this  difference  In  eddy  viscosity  between  the 
two  models  Is  typical  of  the  profiles  within  the 
3-0  Interaction  region.  Within  the  nominal  2-0 
portion  of  the  boundary  layer  upstream  of  the 
interaction,  the  eddy  viscosity  profiles  are  In 
reasonably  close  agreement. 

It  Is  evident  from  Figs.  19  to  22  and  the 
additional  numerous  profiles  studied  that  the 
details  (l.e.,  the  velocity,  pressure  and 
temperature)  of  this  3-0  turbulent  Interaction  are 
relatively  Insensitive  to  the  particular  turbulence 
model  employed,  with  the  exception  of  a  small 
fraction  of  the  boundary  layer  adjacent  to  the 
surface  where  modest  differences  are  observed  In 
the  yaw  angle.  This  implies,  therefore,  that  the 
principal  elements  of  the  flowfleld  structure  are 
rotational  and  Invlscld,  except  near  the  wall  as 
mentioned.  This  represents  a  significant  result 
for  3-0  Interactions,  and  Is  notably  different  from 

2- 0  separated  shock-boundary  layer  Interactions 
wherein  the  differences  between  the  computed 
flowflelds  obtained  using  algebraic  and  two- 
equation  turbulence  models  are  significant4®.  There 
Is  no  reason  to  .  expect,  however,  that  the 
Insensitivity  to  turbulence  model  displayed  In  the 

3- 0  sharp  fin  Interaction  will  necessarily  apply  to 
other  3-0  turbulent  Interactions. 

C.  Flowfleld  Structure  of  3-0  Sharp  Fin 

On  the  basis  of  the  close  agreement  between 
the  computed  and  experimental  .data,  the  computed 
solutions  can  be  utilized  to  examine  the  flow 
structure  of  the  3-0  sharp  fin  Interaction.  The 
close  similarity  of  the  computed  velocity,  yaw  and 
pitch  angle  profiles  for  the  Baldwln-Lomax  and 
Jones-Launder  models  Implies  a  close  agreement  In 
predicted  mean  streamlines,  which  was  confirmed 
through  detailed  comparison  of  nunerous  particle 
traces. 

In  Figs.  23  and  24,  the  computed  surface 
streamlines  ("limiting  streamlines")  obtained  using 
the  Baldwln-Lomax  and  Jones-launder  models  are 
shown.  The  lines  of  coalescence  ("separation")  and 
divergence  ("attachment")  are  observed.  These 
specific  features  are  In  general  agreement  with  the 
experimental  kerosene  lampblack  visualization  as 
described  previously,  although  the  computed  lines 
of  coalescence  appear  further  downstream  than  In 
the  experiment.  In  particular,  the  experimental 
line  of  coalescence  at  z/«*  ■  10  (the  approximate 
spanwlse  limit  of  the  experimental  kerosene 
lampblack  visualization)  Is  xs  *  -5. 7i*,  whereas 
the  computed  lines  of  coalescence  using  the 
Baldwln-Lomax  and  Jones-Launder  models  occur  at  xs 
■  -4.5a.  and  -3.5«a,  respectively.  It  Is  noted 
that  the  discrepancy  between  the  experimental  and 
Baldwln-Lomax  result,  1.2«a,  Is  equal  to  1.2Ax.  In 
the  Jones-Launder  model,  the  discrepancy,  l.e.. 


Z.Ztm  Is  5.6dx.  Further  Investigation  Is  needed  to 
understand  the  effects  of  the  turbulence  models  on 
the  calculated  lines  of  coalescence. 

A  series  of  calculated  mean  streamlines  are 
displayed  in  Figs.  25  to  28,  obtained  from  the 
computed  solution  utilizing  the  3a 1 dwi n-Lomax 
model.  Three  views  of  the  streamlines  are  shown, 
corresponding  to  an  observer  looking  a)  towards  the 
surface  from  above,  b)  downstream,  and  c)  towards 
the  fin  from  the  side.  The  vertical  scale  has  been 
enlarged  by  a  factor  of  three  for  the  purposes  of 
clarity.  In  Fig.  25,  a  series  of  twelve  stream¬ 
lines  are  shown.  Six  streamlines  originate  from 
the  surface,  upstream  of  the  Interaction,  at  equal 
spanwlse  Increments  of  <a,  and  serve  to  define  the 
line  of  coalescence.  Six  additional  streamlines 
originate  Inmedlately  above  the  previous  six,  at  a 
height  of  0. 0048a*.  These  latter  streamlines 
clearly  rise  and  cross  the  line  of  separation,  and 
appear  to  concentrate  within  a  core.  In  Fig.  26, 
another  series  of  twelve  streamlines  are  displayed. 
The  first  six  again  represent  limiting  streamlines, 
and  define  the  line  of  coalescence.  The  second  six 
originate  upstream  at  a  height  of  0.52a*.  These 
particles  display  a  clear  rotational  motion  which 
Is  counterclockwise  as  viewed  looking  downstream. 
In  Fig.  27,  the  three  different  sets  of  six 
streamlines  are  shown.  The  streamlines  display  a 
vortical  structure.  Those  particles  originating 
from  a  higher  y  are  swept  beneath  the  particles 
originating  near  the  surface.  In  Fig.  28,  a  final 
series  of  twelve  streamlines  are  shown.  Again,  a 
set  of  six  surface  streamlines  define  the  line  of 
coalescence.  A  second  set,  originating  upstream  at 
y  »  1.1**,  Is  observed  to  rise  In  the  vicinity  of 
the  line  of  separation,  and  then  drop  towards  the 
surface.  Unlike  the  particles  originating  at  lower 
values  Of  y,  these  particles  eventually  move 
approximately  parallel  to  the  fin  surface. 
Additional  extensive  streamline  patterns  are 
consistent  with  the  above  features. 

A  general  mean  flowfleld  pattern,  developed  on 
the  basis  of  the  streamline  patterns,  is  displayed 
In  Fig.  29.  As  suggested  by  Token^,  the  flowfleld 
Structure  of  the  3-0  sharp  fin  at  the  present 
conditions  Is  dominated  by  a  large  vortical 
structure.  The  line  of  coalescence  (separation) 
defines  the  origin  of  the  3-0  separation  surface 
(Surface  No.  1).  The  line  of  divergence 
(attachment)  represents  the  Intersection  of  a 
second  surface  (Surface  No.  2)  with  the  wall.  This 
second  surface  extends  upstream  into  the. 
undisturbed  flow.  The  fluid  contained  between  the 
wall  and  the  second  surface  is  entrained  into  the 
vortical  structure,  while  the  fluid  above  the 
second  surface  flows  towards  the  wall  and 
approximately  parallel  to  the  fin.  Oue  to  the 
resolution  of  the  nunerlcal  grid,  the  detailed  mean 
flow  structure  In  the  Immediate  vicinity  of  the  fin 
leading  edge  could  not  be  examined.  A  second  small 
vortical  structure  was  observed  within  the  fin 
boundary  layer  and  close  to  the  corner.  In 
agreement  with  the  experimental  observations  of 
Kubota  and  Stollery®  at  Mach  2.3  and  Re{  *  5  x 
104.  No  experimental  data  was  available  for  the 
Immediate  vicinity  of  the  corner  for  this 
configuration,  however,  and  consequently  the 
computed  secondary  vortical  structure  cannot  be 
considered  verified. 


V.  Conclusions 


» 


An  experimental  and  theoretical  study  Is 
performed  for  the  3-0  shock  wave  turbulent  boundary 
layer  Interaction  generated  by  a  sharp  fin  at  Mach 
3  for  a  fin  angle  ag  *20  deg  and  Reynolds  number 
fle4  *  9  x  10= .  Tnls  research  effort  extends 
previous  experimental  and  theoretical  investi¬ 
gations  of  the  3-0  sharp  fin  Interaction  to 
stronger  Interactions.  the  experimental  data 
Include  surface  pressure  profiles,  surface 
streamline  patterns,  and  boundary  layer  profiles  of 
pitot  and  yaw  angle.  Two  separate  theoretical 
approaches  or  "models"  were  employed,  both  of  which 
utilize  the  3-0  compressible  Navler-Stokes 
equations.  The  theoretical  approach  of  Knight 
employs  the  algebraic  turbulent  eddy  viscosity 
model  of  Baldwin  and  Lomax,  and  the  theoretical 
model  of  Horstman  employs  the  two-equation 
turbulence  model  of  Jones  and  Launder  coupled  with 
the  wall  function  model  of  Vlegas  and  Rubesln.  The 
principal  conclusions  are: 

1.  The  computed  surface  pressure,  surface  stream¬ 
lines,  pitot  pressure  and  yaw  angle  profiles 
are  observed  in  good  agreement  with  the 
experimental  data,  thereby  confirming  the 
efficacy  of  the  theoretical  approaches  which 
were  previously  validated  for  the  3-0  sharp  fin 
configuration  at  Mach  3  for  smaller  ag  (l.e. , 
weaker  Interactions). 

2.  The  three  dimensional  velocity  fields  computed 

by  both  models  are  In  close  agreement,  as 
Indicated  by  a  detailed  evaluation  of 
x-component  velocity,  pitch  and  yaw  angle 
profiles.  This  result  was  confirmed  by 
comparison  of  fluid  particle  pathllnes  for 
the  theoretical  models,  which  were  found  to 
be  In  close  agreement.  The  turbulent  eddy 
viscosity  profiles,  however,  differ 

significantly  within  the  3-0  Interaction. 
This  result  implies  that  the  overall  structure 
of  this  3-0  sharp  fin  interaction  is  Insensitive 
to  the  turbulence  model,  except  within  a  small 

•portion  of  the  boundary  layer  adjacent  to  the 
surface.  The  principal  elements  of  the  flow- 
field  structure  are  therefore  rotational  and 
Invlscld,  except  near  the  wall  as  Indicated. 

3.  The  calculated  flowflelds  display  a  prominent 

vortical  structure  associated  with  the  shock¬ 
boundary  layer  Interaction  In  agreement  with 
the  flowfleld  models  of  Token,  and  Kubota  and 
Stollery.  A  three  dimensional  surface  of 

separation  emanates  from  the  line  of 

coalescence  (separation),  and  spirals  Into  the 
vortical  center.  A  second  surface,  emanating 
from  upstream.  Intersects  the  wall  at  the  line 
of  divergence  (attachment),  and  defines  the 
extent  of  the  fluid  entrained  Into  the  vortical 
structure. 
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Fig.  1  Physical  Region  for  3-0  Sharp  Fin 
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Fig.  4  Pitot  Pressure  at  Station  4 
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Fig.  2  Location  of  Experimental  Surveys 


Fig.  5  Pitot  Pressure  at  Station  5 
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Fig.  3  Pitot  Pressure  at  Station  3 


Fig.  6  Pitot  Pressure  at  Station  6 
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Fig.  7  Pitot  Pressure  at  Station  10 


Fig.  10  Pitot  Pressure  at  Station  8 
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Fig.  9  Pitot  Pressure  at  Station  11 


Fig.  12  Yaw  Angle  at  Station  4 


Fig.  15  Yaw  Angle  at  Station  10 


Fig.  18  Yaw  Angle  at  Station  8 


Fig.  27  Computed  Mean  Streamlines  Originating 
Upstream  at  y  *  0,  0.0048a  and 
0.52a.  (Vertical  Scale  Enlarged  by 
Factor  of  3) 
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Fig.  29  Moan  Rowfleld  Structure 


